NRL/MR/6615— 93-6948 


AD-A267  959 


Energy  Deposition  Fluctuations 
Induced  by  Ions  in  Microvolumes  and 
Nanovolumes  —  An  Analytic  Approach 
I.  Theory 

Michael  A.  Xapsos  and  Geoffrey  P.  Summers 
Radiation  Effects  Branch 

Condensed  Matter  &  Radiation  Sciences  Division 

Edward  A.  Burke 

Consultant 
Woburn,  MA 


DTIC 


ELFCTE 
AUG  11 1993 


& 

In 


July  22,  1993 


Approved  lor  public  release;  distribution  unlimited. 


93-17586 


te»‘«93u 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


Public  raportir^  burd«n  for  tNt  coKaction  of  information  «*  oatimatod  to  average  1  hour  per  reaponae,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this 
collection  of  information,  including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  end  Reports,  1216  Jefferson 
Davis  Highway,  Suite  1204.  Arlington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget.  Paperwork  Reduction  Project  (0704-0188),  Washington,  DC  20603. 

1.  AGENCY  USE  ONLY  (Leave  Blank ) 

2.  REPORT  DATE 

3.  REPORT  TYPE  AND  DATES  COVERED 

July  22,  1993 

Final 

4.  TITLE  AND  SUBTITLE 

5.  FUNDING  NUMBERS 

Energy  Deposition  Fluctuations  Induced  by  Ions  in 

Microvolumes  and  Nanovolumes  —  An  Analytic  Approach 

I.  Theory 

TA:  RR-0 12-02-44 

WU:  3824 

6.  AUTHORISI 

Michael  A.  Xapsos,  Edward  A.  Burke,*  and  Geoffrey  P.  Summerst 

7.  PERFORMING  ORGANIZATION  NAMEIS)  AND  ADDRESSIES) 

Naval  Research  Laboratory 

Washington,  DC  20375-5320 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

NRL/MR/6615— 93-6948 

9.  SPONSORING/MONITORING  AGENCY  NAMEIS)  AND  ADDRESSIES) 

Office  of  Naval  Research 

Arlington,  VA  22217-5000 

10.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 

1 1 .  SUPPLEMENTARY  NOTES 

Consultant 

Woburn,  MA  01801 

t  Also  at:  Department  of  Physics 

University  of  Maryland  Baltimore  County 

Baltimore,  MD  21228 

1  2a.  DISTRIBUTION/AVAILABILITY  STATEMENT 

1 2b.  DISTRIBUTION  CODE 

Approved  for  public  release;  distribution  unlimited. 

1 3.  ABSTRACT  ( Maximum  200  words) 

A  new  analytic  method  is  developed  for  calculating  fluctuations  in  energy  deposition  induced  by  ions  in  sites  having  dimensions 
as  small  as  one  nanometer.  The  method  considers  two  contributions  to  the  fluctuation  phenomena.  The  first  is  due  to  direct  ion 
interactions  with  the  site.  The  second  is  due  to  secondary  electrons  that  strike  the  site  but  which  are  produced  by  ion  interactions 
in  the  surrounding  medium.  Particular  attention  is  given  to  the  way  the  two  contributions  are  combined  to  describe  the  complete 
stochastic  process.  This  analytic  combination  of  the  two  contributions  has  not  previously  been  developed  in  the  context  of 
microdosimetry  where  Monte  Carlo  simulations  are  often  used  to  describe  the  combined  processes.  The  new  analytic  approach  is 
much  more  versatile  than  Monte  Carlo  methods  and  can  easily  be  applied  to  different  incident  radiations,  target  materials  and  target 
geometries.  In  addition  to  its  versatility  and  convenience,  the  physical  factors  which  contribute  to  the  fluctuations  in  energy 
deposition  are  brought  clearly  into  focus.  Furthermore,  the  method  allows  calculations  of  ionization  fluctuations  to  be  performed. 

The  purpose  of  this  report  is  to  present  detailed  derivations  of  the  basic  equations  of  the  model  beginning  with  concepts  that  are 
well  known  in  stochastic  theory.  Additionally,  general  expressions  for  the  important  microdosimetric  quantities  associated  with  the 
frequency  and  dose  weighted  averages  of  specific  energy  and  lineal  energy  are  obtained.  These  microdosimetric  quantities  are 
useful  for  analyzing  energy  distributions  and  event  frequencies  in  small  sites.  Emphasis  is  given  to  explaining  concepts  used  almost 
exclusively  in  the  field  of  microdosimetry  to  an  unfamiliar  reader.  A  later  report  will  show  how  the  model  is  applied  and  the 
results  that  are  obtained  in  actual  examples. 

14.  SUBJECT  TERMS 

15.  NUMBER  OF  PAGES 

Energy  deposition  fluctuations 

Nanovolumes 

66 

Microdosimetry 

Microvolumes 

16.  PRICE  CODE 

1 7.  SECURITY  CLASSIFICATION 

OF  REPORT 

16.  SECURITY  CLASSIFICATION 

OF  THIS  PAGE 

19.  SECURITY  CLASSIFICATION 

OF  ABSTRACT 

20.  LIMITATION  OF  ABSTRACT 

UNCLASSIFIED 

UNCLASSIFIED 

UNCLASSIFIED 

UL 

NSN  7640-01-280-5500 

i 

Standard  Form  298  Iftav.  2-89) 
Praacribad  by  ANSI  Std  239-18 

298-102 

CONTENTS 


% 

NOTATION  .  v 

I.  INTRODUCTION .  1 

II.  THEORY  .  2 

A.  Review  of  the  Statistics  of  Continuous  Random  Variables .  2 

B.  Relations  between  the  Random  Physical  Variables  and 

Energy  Deposition  Fluctuations .  4 

1.  Random  Physical  Variables  in  the  Energy  Deposition  Process  .  5 

2.  Energy  Deposition  Fluctuations  with  Random  Variables 

Described  by  Moments . 10 

a.  Dependence  of  Dose  Fluctuations  on  Average  Dose  . 10 

b.  Fluctuations  in  Dose  due  to  Single  Ion  Events . 16 

3.  Energy  Deposition  Fluctuations  with  Random  Variables  Described 

by  Relative  Variances . 21 

a.  Relations  between  Moment  Ratios  and  Relative  Variances  of 

the  Random  Variables  . 21 

b.  Conversion  to  the  Relative  Variance  Description  . 23 

C.  The  Two  Component  Model . 26 

1.  Combined  or  Mixed  Probability  Distributions  . 26 

2.  Energy  Deposition  Fluctuations  with  Random  Variables 

»  Described  by  Moments . 30 

3.  Energy  Deposition  Fluctuations  with  Random  Variables  Described  by 

*  Relative  Variances . 37 


m 


D.  Microdosimetric  Quantities  and  the  Two  Component  Model  . 41 

1.  Specific  Energy  Related  Quantities . 41 

a.  Definition  and  Background . 41 

b.  Specific  Energy  and  the  Two  Component  Model  . 44 

2.  Lineal  Energy  Related  Quantities . 46 

a.  Definition  and  Background . 46 

b.  Lineal  Energy  and  the  Two  Component  Model  . 49 

E.  Overview . 53 

ID.  SUMMARY . 57 

REFERENCES  . 58 


Accesion  For  ( 

NTIS 

CRA-Sil 

td 

DTIC 

TAB 

Unan;) 

JjStliK 

By 

|  Di-sMb 

-  t'O::  / 

Av_i, :aO'T. 

•/  Codes 

Avan 

*  i  /  or 

Dist 

S\.i. 

cijl 

0 

DTIC  QUALITY  INSPECTED  3 


IV 


NOTATION 


F  average  chord  length  of  the  site 
D  total  dose  deposited  in  the  site 
De  dose  deposited  in  the  site  by  electrons 

Dion  dose  deposited  in  the  site  by  ions 

d  dose  deposited  in  the  site  by  one  particle 

de  dose  deposited  in  the  site  by  one  electron 

dion  dose  deposited  in  the  site  by  one  ion 
d(y)  dose  weighted  probability  density  function  for  lineal  energy 
d(z)  dose  weighted  probability  density  function  for  specific  energy 

d{(z)  dose  weighted  probability  density  function  for  specific  energy  per  event 

E  total  energy  deposited  in  the  site 
Ee  energy  deposited  in  the  site  by  electrons 

Eion  energy  deposited  in  the  site  by  ions 

F  Fano  factor 

Fe  Fano  factor  for  electrons 

Fion  Fano  factor  for  ions 

fon  fraction  of  dose  deposited  in  the  site  by  ions 

1  -fion  fraction  of  dose  deposited  in  the  site  by  electrons 

f(N)  probability  density  function  for  number  of  particle  events  in  the  site 
f(x)  probability  density  function  for  number  of  occurrences  of  x 
f(y)  frequency  weighted  probability  density  function  for  lineal  energy 
ffz)  frequency  weighted  probability  density  function  for  specific  energy 

/,  (z)  frequency  weighted  probability  density  function  for  specific  energy  per  event 

L  (unrestricted)  Linear  Energy  Transfer  of  a  particle  in  the  site 
Le  (unrestricted)  Linear  Energy  Transfer  of  an  electron  in  the  site 

Lion  (unrestricted)  Linear  Energy  Transfer  of  an  ion  in  the  site 

N  number  of  events  in  the  site 
N'  fixed  value  of  N 

Na  number  of  events  in  the  /1-distribution 

Nb  number  of  events  in  the  fi-distribution 

Ne  number  of  electron  events  in  the  site 
Nion  number  of  ion  events  in  the  site 
n  number  of  collisions  undergone  by  a  particle  in  the  site 
ne  number  of  collisions  undergone  by  an  electron  in  the  site 

nion  number  of  collisions  undergone  by  an  ion  in  the  site 

P  normalized  or  relative  frequency  of  the  /1-distribution;  fraction  of  events  due  to  ions 
\-P  normalized  or  relative  frequency  of  the  fi-distribution;  fraction  of  events  due  to  electrons 
S0  surface  area  of  the  site 

s  path  length  of  a  particle  in  the  site 
se  path  length  of  an  electron  in  the  site 
sion  path  length  of  an  ion  in  the  site 
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VD  relative  variance  of  the  total  dose  deposited  in  the  site 

VD  e  relative  variance  of  the  dose  deposited  in  the  site  by  electrons 

VD  ion  relative  variance  of  the  dose  deposited  in  the  site  by  ions 

VE  relative  variance  of  the  total  energy  deposited  in  the  site 

VF  relative  variance  of  Fano  fluctuations 

VFe  relative  variance  of  Fano  fluctuations  due  to  electrons 

VF  ion  relative  variance  of  Fano  fluctuations  due  to  ions 

VL  relative  variance  of  the  LET  distribution  in  the  site 

VL  e  relative  variance  of  the  electron  LET  distribution  in  the  site 

V L  ion  relative  variance  of  the  ion  LET  distribution  in  the  site 

VN  relative  variance  of  the  number  of  events  in  the  site 

Vs  relative  variance  of  the  path  length  distribution  in  the  site 

Vse  relative  variance  of  the  electron  path  length  distribution  in  the  site 

Vs  ion  relative  variance  of  the  ion  path  length  distribution  in  the  site 

Vs!r  relative  variance  of  the  energy-loss  straggling  distribution  in  the  site 

Vstre  relative  variance  of  the  energy-loss  straggling  distribution  of  electrons  in  the  site 

VSIr  jon  relative  variance  of  the  energy-loss  straggling  distribution  of  ions  in  the  site 

V6  relative  variance  of  the  energy  deposited  in  the  site  by  one  collision 

V£  relative  variance  of  the  energy  deposited  in  the  site  by  one  particle 

V0  volume  of  the  site 

V,  single  event  relative  variance 

Vl  e  single  event  relative  variance  of  electrons 

Vl  lon  single  event  relative  variance  of  ions 

varF(e)  variance  of  Fano  fluctuations 

varL.s(e)  variance  of  e  due  to  L  and  s 

varslr(e)  variance  of  e  due  to  energy-loss  straggling 

vv  average  energy  required  to  produce  one  electron-hole  pair  in  the  site 

We  average  energy  required  for  electrons  to  produce  one  electron-hole  pair  in  the  site 

vv“  average  energy  required  for  ions  to  produce  one  electron-hole  pair  in  the  site 

w(x)  probability  density  function  for  a  weighted  distribution 

x  general,  continuous  random  variable 

xA  continuous  random  variable  of  the  /4-distribution 

xB  continuous  random  variable  of  the  ^-distribution 

x^  weighted  average  of  x 

y  lineal  energy 

yP  dose  mean  lineal  energy 

y^e  dose  mean  lineal  energy  for  electrons 

Jo'jan  dose  mean  lineal  energy  for  ions 

yp  frequency  mean  lineal  energy 

yp~  frequency  mean  lineal  energy  for  electrons 

yp~~n  frequency  mean  lineal  energy  for  ions 

z  specific  energy  (in  section  Dla);  otherwise,  specific  energy  per  event 
dose  mean  specific  energy  per  event 


vi 


~  dose  mean  specific  energy  per  electron  event 
dose  mean  specific  energy  per  ion  event 
zf  frequency  mean  specific  energy  per  event 
z^~  frequency  mean  specific  energy  per  electron  event 
zFJOn  frequency  mean  specific  energy  per  ion  event 
6  energy  deposited  by  a  particle  in  the  site  by  one  collision 
6e  energy  deposited  by  an  electron  in  the  site  by  one  collision 
6ion  energy  deposited  by  an  ion  in  the  site  by  one  collision 

e  energy  deposited  in  the  site  by  one  particle 

te  energy  deposited  in  the  site  by  one  electron 

tion  energy  deposited  in  the  site  by  one  ion 

ti  number  of  electron-hole  pairs  produced  in  the  site  by  one  particle 

r\e  number  of  electron-hole  pairs  produced  in  the  site  by  one  electron 

x\ion  number  of  electron-hole  pairs  produced  in  the  site  by  one  ion 

p0  density  of  the  material  composing  the  site 
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ENERGY  DEPOSITION  FLUCTUATIONS  INDUCED  BY  IONS  IN 
MICRO  VOLUMES  AND  NANO  VOLUMES  -  AN  ANALYTIC  APPROACH 

I.  THEORY 

I.  INTRODUCTION 

The  spatial  distribution  of  radiation  energy  deposited  on  a  micrometer  to  nanometer  scale 
is  central  to  understanding  radiation  effects  in  biological  systems1,2  and  modem  microelectronics.3 
Because  of  the  essential  granular  nature  of  all  radiation,  energy  deposition  in  small  volumes  is 
a  statistical  process  so  that  knowledge  of  the  macroscopically  measured  dose  does  not  determine 
the  energy  imparted  to  a  particular  small  volume.  Monte  Carlo  simulations  have  been  the 
traditional  method  to  model  energy  deposition  by  ions  in  sites  whose  dimensions  are  less  than 
approximately  one  micrometer.  Although  such  simulations  have  proven  to  be  valuable,  they 
involve  sophisticated  computer  codes  and  time  consuming  scoring  techniques,  require  large 
quantities  of  input  information,  tend  to  be  inflexible,  and  are  often  tedious  to  interpret  in  terms 
of  the  relative  effects  of  various  input  parameters.  Consequently,  certain  aspects  of 
radiobiological  and  microelectronic  studies  are  difficult  to  address  with  Monte  Carlo  calculations. 
On  the  other  hand,  an  analytical  approach  can  provide  information  to  supplement  Monte  Carlo 
results,  and  is  more  suited  to  analyze  the  wide  ranges  of  conditions  often  encountered  in  radiation 
studies. 

Experiments  in  radiobiological  research,  for  example,  are  commonly  performed  with  a 
variety  of  ions  and  energies  to  cover  a  broad  range  of  Linear  Energy  Transfer  (LET).  There  are 
also  uncertainties  associated  with  the  dimensions  and  geometry  of  the  sensitive  site.  The 
dimensions  may  vary  anywhere  from  the  size  of  a  cell  nucleus  (micrometers)  to  that  of  a  DNA 
molecule  (nanometers),  and  the  geometry  is  determined  accordingly.  By  contrast,  Monte  Carlo 
calculations  have  been  mainly  limited  to  proton  irradiations  of  sites  with  a  spherical  geometry. 
Furthermore,  in  specific  cases  such  as  radiation-induced  DNA  damage,  it  is  not  generally  agreed 
upon  whether  energy  deposition  (ionization  plus  excitation)  or  simply  electron-hole  pair 
production  (ionization)  is  the  relevant  quantity  related  to  the  damage  mechanism.  Adjustments 
of  this  nature  in  Monte  Carlo  codes  are  difficult  to  make. 

The  situations  which  occur  in  radiation  effects  research  on  microelectronics  have  related 
concerns.  A  variety  of  ions  are  commonly  used  in  experiments,  and  the  sensitive  volumes  of 
targets  involve  different  materials  (Si,  GaAs,  SiO^,  dimensions  (nanometers  to  10’s  of 
micrometers)  and  geometries  associated  with  depletion,  diffusion  and  insulating  regions  of 
devices.  It  is  clear  that  there  is  a  need  for  a  general  method  which  can  quickly  determine  the 
effects  that  basic  input  parameters  have  on  the  overall  outcome  of  irradiating  very  small  volumes 
of  matter  with  ions.  This  work  presents  such  a  method  for  calculating  energy  deposition 
fluctuations  which  is  analytical,  based  on  first  principles  and  whose  basic  equations  do  not 
involve  fitted  parameters. 

The  basic  derivations  are  presented  in  Section  II.  Some  background  in  statistical 
terminology  is  reviewed  in  IIA,  and  then  a  new  derivation  of  relations  introduced  by  Kellerer 
is  presented  in  IIB.  This  relates  the  energy  deposition  fluctuations  caused  by  a  single  type  of 
event  to  the  first  and  second  moments  of  three  fundamental  distributions.  In  IIC,  the  relations 
are  generalized  to  the  "Two  Component  Model",  which  analytically  describes  ion  induced  energy 
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deposition  fluctuations  for  dimensions  as  small  as  one  nanometer.  In  these  sections,  the  methods 
are  presented  in  terms  of  standard  statistical  quantities  such  as  variances,  relative  variances  and 
moments  of  distributions.  The  microdosimetry  literature,  however,  has  evolved  using  somewhat 
different  terminology.  In  IID,  the  relation  of  this  terminology  to  discussion  in  previous  se<  ’ons 
is  shown.  Those  readers  not  interested  in  this  section  can  simply  omit  it  without  any  I  ss  of 
continuity.  Section  IIE  gives,  with  the  aid  of  flow  charts,  an  overview  of  the  ideas  behind  the 
Two  Component  Model.  Finally,  section  IH  is  a  brief  summary. 


H.  THEORY 

A.  Review  of  the  Statistics  of  Continuous  Random  Variables 

In  much  of  what  follows,  it  will  be  useful  to  parameterize  the  distribution  of  values  which 
a  random  variable  can  assume  in  terms  of  its  relative  variance  and  in  terms  of  the  moments  of 
the  distribution.  Suppose  x  is  a  continuous  random  variable  described  by  a  probability  density 
function  f(x)  that  gives  the  probability  of  the  number  of  occurrences  of  x.  Then  f(x)dx  represents 
the  probability  that  the  random  variable  x  is  between  x  and  x  +  dx.  The  first  moment  of  x,  the 
average,  is  given  by 


x  =  fx‘f(x)dx  (D 

where  the  overlining  of  x  indicates  it  is  an  average.  In  general,  the  m-  moment  of  x  is  given  by 

•• 

x m  =  Jxm  'fix)dx  (2) 


The  variance  of  x,  denoted  by  var(x),  is  given  by 


var(x)  =  o\  =  f(x  -  x)2-fc)dx 


(3) 


where  ax  is  the  standard  deviation  of  x.  The  variance  is  also  commonly  written  as  the  second 
moment  minus  the  square  of  the  first  moment. 
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(4) 


var{x)  =  x1  -  (x)2 


The  relative  variance  of  the  random  variable  x ,  Vx,  is  defined  as  the  variance  divided  by  the 
square  of  the  average. 


var(x) 

(x)2 


(5) 


Another  quantity  commonly  used  to  characterize  the  variability  of  x  is  the  coefficient  of 
variation,  ajx,  which  is  the  square  root  of  the  relative  variance. 

Rather  than  use  the  average  and  variance  (or  relative  variance)  to  describe  the  distribution 
of  values  that  x  can  assume,  it  is  sometimes  more  convenient  to  use  the  first  and  second  moments 
of  the  distribution.  A  useftil  relation  to  know  is  that  between  the  ratio  of  the  second  moment  to 
the  first  moment  of  x  and  the  relative  variance  of  x.  Dividing  each  side  of  equation  (4)  by  the 
square  of  the  average  of  x,  and  making  use  of  equation  (5),  the  moment  ratio  is  easily  shown  to 
be  given  by 


x 


x 


•(i  ♦  vx) 


(« 


There  are  many  other  types  of  probability  density  functions  besides  f(x)  that  can  be 
employed.  One  that  will  be  useful,  denoted  here  as  w(x),  describes  the  probability  that  the 
variable  x  is  a  certain  value,  i.e., 


w(x)dx 


x-Mdx 


fx-ftx)dx 


(7) 


In  other  words,  w(x)dx  is  the  fraction  of  x  associated  with  a  particular  value  of  x.  This  type  of 
distribution  is  referred  to  as  a  weighted  distribution.  It  should  be  contrasted  to  the  distribution 
given  by  f(x)  which  describes  the  number  of  occurrences  of  x.  The  weighted  distribution,  w(x), 
is  useful  for  analyzing  the  concentration  of  x  rather  than  its  frequency  of  occurrence. 

The  average  of  x  taken  over  the  weighted  distribution  is 
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xw  =  jx  *w(x )dx 


(8) 


where  the  subscript  W  is  used  to  indicate  it  is  a  weighted  average  to  distinguish  it  from  the 
arithmetic  average  given  by  equation  (1).  The  w  hted  average  can  be  expressed  in  terms  of 
f(x )  by  substituting  equation  (7)  into  equation  (8). 


/  x2-fa)dx 


f  x  'fix)dx 


(9) 


It  will  be  recognized  from  equation  (2)  that  this  weighted  average  can  be  expressed  as  the  second 
moment  of  x  divided  by  the  first  moment  of  x. 


x 


w  ~ 


(10) 


Expressions  involving  the  ratio  of  the  second  moment  to  the  first  moment  of  a  continuous 
random  variable  play  a  fundamental  role  in  microdosimetry  theory  because  the  concentration  of 
the  physical  quantity  which  the  variable  represents  can  influence  the  outcome  of  an  irradiation. 


B.  Relations  between  the  Random  Physical  Variables  and  Energy  Deposition  Fluctuations 

The  basic  relations  describing  energy  deposition  fluctuations  for  a  single  type  of  event 
such  as  the  traversal  of  a  target  site  by  an  ion  was  introduced  originally  by  Kellerer.1  However, 
the  origin  of  the  equations  is  only  discussed  briefly  in  this  somewhat  obscure  reference.  One 
problem  confronting  a  reader  not  familiar  with  the  theory  is  its  apparent  isolation  from  other 
fields  of  research.  The  reason  for  presenting  the  following  derivation  is  to  make  clear  the 
principles  behind  the  equations  so  that  the  model  which  is  developed  from  these  relations  can  be 
applied  to  other  physical  systems.  To  this  end,  a  more  straight-forward  derivation  of  the  basic 
equations  than  that  outlined  in  reference  1  is  presented. 

In  section  Bl,  a  qualitative  description  of  the  random  variables  important  in  the  energy 
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deposition  process  is  presented.  In  B2,  the  results  introduced  by  Kellerer  are  obtained  from 
concepts  which  are  well  known  in  stochastic  theory.  These  equations  have  the  random  variables 
expressed  in  terms  of  the  moments  of  their  distributions.  In  section  B3,  the  approach  is 
converted  to  an  equivalent  one  in  which  the  random  variables  are  described  by  the  relative 
variance  of  the  appropriate  distribution. 


1.  Random  Physical  Variables  in  the  Energy  Deposition  Process 

Figure  la  gives  an  overall  view  of  the  energy  deposition  process.  Ions  are  considered 
to  traverse  a  small  volume  of  material  in  straight  paths  without  losing  a  significant  portion  of 
their  energy  in  the  volume.  The  paths  are  indicated  by  arrows  and  the  particular  situation  shown 
is  for  a  site  with  an  arbitrary  geometry  crossed  by  three  particles.  More  generally,  the  number 
of  particles  which  traverse  the  site,  N,  is  subject  to  statistical  fluctuations.  N  is  also  referred  to 
as  the  number  of  events  which  occur  in  the  site.  The  particles  are  assumed  to  arrive  at  the  site 
independently  ar:d  the  variability  of  N  is  therefore  described  by  Poisson  statistics.  More 
significantly,  the  energy  imparted  to  the  volume  by  a  single  particle,  e,  also  displays  variability 
and  will  be  discussed  below.  The  overall  process  is  therefore  compound  Poisson,  i.e. ,  a  process 
of  N  independent  random  events  of  varying  magnitude  e. 

Next  the  random  physical  variables  associated  with  a  single  particle  traversal  of  the  site 
are  considered.  As  figure  lb  indicates,  the  particles  can  have  a  distribution  of  LETs  as  well  as 
a  distribution  of  path  lengths  through  the  site.  The  former  is  determined  by  the  type  and  energy 
of  the  incident  particle  while  the  latter  depends  on  the  geometry  of  the  site  and  the  directions 
from  which  the  particles  are  incident.  The  energy  deposited  by  a  particle  traversing  the  site 
depends  on  the  product  of  the  random  variables  L  and  s.  L  is  the  LET  of  the  particle  in  the 
material  under  consideration  and  s  is  the  particle’s  path  length  through  the  site. 

The  next  process  to  consider  results  from  fluctuations  in  the  energy  deposition  for  a  single 
particle  traversing  the  site  with  a  fixed  LET  and  a  f  xed  path  length.  This  process  is  referred 
to  as  energy-loss  straggling.  As  the  particle  crosses  the  site,  it  undergoes  a  random  number  of 
collisions  with  the  electrons  in  the  site.  The  number  of  collisions  or  energy  transfer  points 
depicted  in  figure  lc  is  four  and  in  general  is  denoted  by  n.  The  collisions  occur  independently, 
so  n  is  described  by  Poisson  statistics.  Each  collision  has  a  magnitude  described  by  the  random 
variable  6,  which  accounts  for  both  ionizations  and  inelastic  excitations.  The  energy-loss 
straggling  is  also  a  compound  Poisson  process  because  it  consists  of  n  independent  events  of 
varying  magnitude  6.  The  compound  Poisson  process  therefore  plays  a  role  at  two  levels  of  the 
energy  deposition  process,  as  has  been  pointed  out  previously  by  Kellerer.2 

The  above  discussion  considers  the  random  factors  contributing  to  the  energy  deposited 
in  a  small  volume  of  matter  which  leads  in  general  to  both  ionizations  and  inelastic  excitations. 
This  may  be  the  quantity  of  interest,  for  example,  in  considering  radiation-induced  damage  in 
DNA  or  other  polymers.  However,  a  distinction  must  be  made  from  the  situation  where  the 
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Figure  la.  Energy  deposition  in  a  small  site. 
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L  RANDOM 
S  RANDOM 

Figure  lb.  LET  and  path  length  random  variables. 
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Figure  lc.  Energy-loss  straggling. 
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Figure  Id.  Fano  fluctuations. 
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relevant  process  involves  only  ionization.  This  latter  situation  does  not  depend  on  inelastic 
excitations  and  is  relevant  to  many  types  of  measurements.  Examples  include  proportional 
counter  spectra  used  in  the  microdosimetry  field  and  current-voltage  measurements  of  MOS 
transistors  at  low  temperatures,  which  yield  information  about  the  distribution  of  radiation- 
generated  holes  in  the  oxide  layer. 

Figure  Id  shows  the  random  factors  involved  in  converting  energy  deposited  in  the 
volume  to  ionizations.  Consider  an  ion  crossing  the  site  and  depositing  energy  e.  Each 
interaction  in  the  site  may  result  in  electron-hole  pair  production  (indicated  by  the  paired  "  +  " 
and  signs),  inelastic  excitations  (denoted  by  asterisks,*),  or  both.  The  statistical  partitioning 
of  the  energy  results  in  a  random  number  of  electron-hole  pairs  being  produced  in  the  site,  ti  . 
This  variability  is  referred  to  as  Fano  fluctuations.4  The  electron-hole  pairs  are  not  produced 
independently,  so  the  process  is  not  Poisson.  Fano  fluctuations  will  be  considered  in  more  detail 
in  section  B2b. 

To  summarize  this  section,  the  discussion  involving  figure  la  gives  an  overview  of  the 
entire  process.  Figures  lb  and  lc  indicate  the  variability  of  the  energy  deposited  by  a  single 
particle.  Finally,  figure  Id  indicates  the  additional  variability  that  must  be  considered  to  convert 
energy  deposited  by  a  single  particle  to  ionization. 


2,  Energy  Deposition  Fluctuations  with  Random  Variables  Described  by  Moments 

The  overall  goal  of  section  B2  is  to  derive  the  relations  introduced  by  Kellerer  which 
describe  the  energy  deposition  fluctuations  in  a  site  traversed  by  an  incident  particle.  For  the 
situation  when  the  incident  particle  is  an  ion,  the  events  are  commonly  referred  to  in  the 
microdosimetry  literature  as  ion  events,  direct  events  or  "crossers".  The  equations  depend  on 
the  random  variables  discussed  in  section  Bl.  The  first  basic  equation,  to  be  derived  in  B2a, 
describes  the  dependence  of  the  fluctuations  on  the  average  (or  macroscopically  measured)  dose 
received  by  the  medium.  The  second  basic  equation,  obtained  in  B2b,  describes  the  fluctuations 
in  the  energy  deposited  by  a  single  particle  traversing  the  site. 


a.  Dependence  of  Dose  Fluctuations  on  Average  Dose 

Assume  the  target  site  is  struck  by  a  random  number  of  incident  ions,  N.  The  average 
(or  expectation  value)  of  N  is 


N  =  YlN-m 

AM) 


(ID 


where  f(N)  is  the  probability  that  N  ions  cross  the  site.  Note  that  f(N)  is  analogous  to  the 
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probability  density  function  in  equation  (1)  except  that  here  N  is  a  discrete  random  variable. 


The  total  energy  deposited  in  the  site,  E,  is  a  random  variable  given  by 

E  =  Cj  +  e2  +  -  +  tff  (12) 


where  e,  denotes  the  energy  deposited  by  the  /-  of  N  total  ions.  It  is  assumed  that  each  ion 
arrives  independently.  The  e,  are  therefore  mutually  independent,  identically  distributed  random 
variables.  From  equation  (4),  the  variance  of  the  energy  deposited  by  a  single  particle  traversing 
the  site  is 


var(e)  *  e2  -  (e)2 


(13) 


The  initial  objective  in  this  section  is  to  determine  the  average  energy  deposited  in  the  site  and 
its  variance  in  terms  of  the  random  variables  e(  and  N. 


First,  an  expression  for  the  average  energy  deposited  will  be  obtained.  Since  this  quantity 
depends  on  both  zt  and  N,  it  will  be  determined  by  first  fixing  N  at  an  arbitrary  value  of  rf.  The 
average  energy  is  said  to  be  conditional  on  N.  Its  conditional  expectation  is  written 


(14) 


The  left-hand  side  of  equation  (14)  denotes  the  average  value  of  E  with  N  held  at  rf.  Its  value 
is  obtained  by  summing  the  average  energy  deposited  in  the  site  by  particles  1  through  IV'.  The 
random  variables  et  are  independent  and  identically  distributed  so, 


e 


i*  l  = 


ci 


e 


(15) 


That  is  all  the  z,  are  equal.  Equation  (14)  can  then  be  written 
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=  ^/,£ 


(16) 


Now  the  expression  for  the  average  total  energy,  E,  can  be  obtained  by  weighting  its  conditional 
expectation,  and  summing  over  all  possible  values  of  N'.  This  summation  gives 


e  =  £ N'- £  -w') 

N1- 0 


(17) 


Since  e  is  independent  of  rf,  it  can  be  taken  outside  the  summation.  Furthermore,  A/  is  simply 
a  dummy  variable  for  N.  Making  use  of  these  facts  and  equation  (11)  gives  the  result 

E  =  t-N  (18) 


This  indicates  the  expected  result  that  the  average  energy  deposited  in  the  site  is  the  product  of 
the  average  energy  deposited  per  event  and  the  average  number  of  events. 

Next  the  variance  of  the  energy  deposited  in  the  volume  will  be  determined.  With  the 
average  value  of  E  now  known,  its  variance  can  be  obtained  using 

W  =  var(E)  +  (E)2  (19) 


once  the  second  moment  of  E  is  determined.  Using  an  analogous  procedure  to  that  above,  N  is 
set  equal  to  a  fixed  A/  to  obtain  the  conditional  expectation  of  the  average  of  the  square  of  E. 

e\.N'  =  +  (e)\.n>  (20) 

The  second  term  on  the  right  hand  side  of  equation  (20)  is  known  from  (16).  The  first  term  can 
be  evaluated  making  use  of  the  properties  that  the  e,  are  independent  and  identically  distributed 
random  variables.  First  of  all,  because  the  e,  are  independent,  the  variance  contributed  by  INK 
particles  is  the  sum  of  the  individual  variances,  i.e., 
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(21) 


N' 


v«K£)U.y  ■  EwK«i) 

t-i 


Furthermore,  because  the  e,  are  identically  distributed,  (21)  simplifies  to 


vor(£) \NmN>  =  N'- vatic) 


(22) 


Substituting  the  results  of  equations  (16)  and  (22)  into  the  right-hand  side  of  (20), 


E\.n>  =  N' vatic)  +  (Nfic)2 


(23) 


As  before,  the  conditional  expectation  must  be  weighted  by  /(TV)  and  summed  over  all  possible 
values  of  rf  to  obtain  the  average.  The  average  of  the  square  of  the  energy  deposited  in  the  site 
is  therefore  given  by 


£2  =  £  [tf'-voKe)  +  {Nrf‘(t)2]‘j(Nr)  (24) 

N'-O 


Equ^on  (24)  is  easily  simplified  to  give 

E2  =  vatic) N  +  (cf-N2 


(25) 


Substituting  the  expressions  for  the  first  and  second  moments  of  E  given  by  (18)  and  (25), 
respectively,  into  (19)  gives  the  variance  of  E. 


var(E)  =  vatic)  ‘  N  +  (c)2*^  -  (IV)2] 


(26) 


From  the  definition  of  variance  given  by  equation  (4),  equation  (26)  can  also  be  written 
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var(E)  =  var(z)‘N  +  ( z)2-var(N ) 


(27) 


This  result  is  equivalent  to  Burgess’  variance  theorem  or  simply  the  variance  theorem5,  which 
is  well  known  in  stochastic  theory  as  the  basic  expression  for  determining  the  variance  of  a  sum 
of  a  random  number  of  random  variables.* 

As  discussed  in  section  Bl,  the  number  N of  random  variables  e  which  make  up  the  sum 
in  equation  (12)  is  controlled  by  Poisson  statistics.  A  variable  controlled  by  Poisson  statistics 
has  die  property  that  its  variance  is  equal  to  its  average  value,  so 


var(N)  =  N 


(28) 


Substituting  (28)  into  (27)  results  in 


var{E)  =  N  \var{z)  +  (e)2] 


(29) 


Using  equation  (13),  this  can  be  further  simplified  to  give 


var(E)  «  N-z2 


(30) 


The  initial  objective  is  now  complete.  Equations  (18)  and  (30)  give  the  average  and  the  variance 
of  the  energy  deposited  in  the  site,  respectively,  in  terms  of  the  random  variables  N  and  c.  It 
should  be  noted  that  the  particularly  simple  form  of  (30)  is  a  direct  result  of  the  fact  that  energy 
deposition  in  a  small  site  is  a  compound  Poisson  process.  The  total  variance  of  such  a  process 
depends  only  on  the  average  number  of  events  and  the  second  moment  of  the  distribution  of 
single  events. 

Next,  these  results  will  be  used  to  determine  the  variability  of  the  dose  deposited  in  the 
small  volume  as  a  function  of  the  average  dose.  Dividing  equation  (30)  by  die  square  of 
equation  (18)  gives  the  relative  variance  of  E,  VE. 


*  Another  equivalent  form  of  Burgess’  variance  theorem  can  be  obtained  by  converting  the 
three  variances  in  equation  (27)  to  relative  variances.  The  result  is  VE  =  VN  +  VJN. 
Breitenberger  (reference  1 1)  gives  this  form  of  the  result  for  what  he  refers  to  as  two  processes 
in  "cascade". 
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(31) 


y  =  varjE)  =  e2 
£  ~  (E)2  (t)2  N 


Using  (18)  to  eliminate  the  dependence  on  N  results  in 


(32) 


In  evaluating  (32),  e  and  E  must  have  the  same  energy  units.  In  microdosimetry 
calculations,  e  typically  is  in  keV.  However,  the  energy  deposited  in  the  site,  E,  is  commonly 
referred  to  in  terms  of  radiation  absorbed  dose  (rad)  units.  (1  rad  =  .01  Gray  =  100  erg/g) 
The  appropriate  conversion  is: 


=  1.6xl0~u 


(33) 


where  E  is  in  keV  and  D  is  the  corresponding  dose  in  rads.  V0  is  the  volume  of  the  site  in  cm3 
and  p0  is  the  density  in  g/cm3.  Since  D  is  proportional  to  E,  it  follows  that  the  relative  variance 
of  D,  VD,  equals  the  relative  variance  of  E,  VE.  Using  this  equivalence  in  equation  (32),  and 
making  use  of  equation  (33)  to  express  E  in  terms  of  rad  units,  equation  (32)  becomes 


1  . 1.6xl0~n  ,  e2 
D  *Vp0  e 


(34) 


where  e  is  in  keV  and  D  in  rads.  VD  is  the  relative  variance  of  the  dose  deposited  in  the  site. 
Furthermore,  recall  that  c  is  the  energy  deposited  in  the  site  by  a  single  ion,  which  is 
independent  of  dose.  The  right  hand  side  of  (34)  is  therefore  the  product  of  a  dose  independent 
factor  and  the  reciprocal  of  the  average  dose.  Denoting  the  dose  independent  quantity  by  this 
can  be  written 


(35) 
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where  z£  is  given  by  (36)  below.  Equation  (35)  is  the  final  result  of  this  section.  It  describes 
the  relative  variance  of  the  dose  deposited  in  a  small  volume  as  a  function  of  the  average  or 
macroscopically  measured  dose.  As  noted  in  IIA,  a  relative  variance  can  be  converted  to  a 
variance  by  multiplying  by  the  square  of  the  average  or  to  a  coefficient  of  variation  by  taking 
its  square  root.  The  quantity  z p  will  be  discussed  in  detail  in  the  following  section. 

The  main  results  of  this  section  can  be  summarized  as  follows.  Equations  (18)  and  (27) 
give  the  average  and  variance  of  E,  respectively,  in  terms  of  the  random  variables  e  and  N. 
Equation  (27)  is  also  equivalent  to  Burgess’  variance  theorem  or  the  variance  theorem.  If  it  is 
assumed  that  the  variable  N  is  controlled  by  a  Poisson  process,  then  the  more  general  equation 
(27)  reduces  to  equation  (30),  which  describes  a  compound  Poisson  process.  Equation  (30)  can 
be  algebraically  manipulated  to  obtain  equation  (35),  which  expresses  the  relative  variance  of  the 
dose  deposited  in  a  small  site  in  terms  of  the  average  dose. 


b.  Fluctuations  in  Dose  due  to  Single  Ion  Events 

Comparing  equations  (34)  and  (35)  it  is  seen  that  the  dose  independent  quantity  denoted 
by  zZ  is  given  by 


zd  ~ 


1.6x10'“ .  e2 
Vo  Po  e 


(36) 


Equation  (36)  can  also  be  written  as 


= 


£ 

d 


(37) 


where  d  is  the  dose  deposited  in  the  site  by  a  single  ion.  In  the  microdosimetry  literature,  TD  is 
called  the  dose  mean  specific  energy  per  event.6  It  is  the  ratio  of  the  second  moment  to  the  first 
moment  of  the  dose  deposited  in  the  site  due  to  single  events.  Note  from  equation  (10)  that  it 
can  be  regarded  as  a  weighted  average  of  the  dose  deposited  by  single  events.  As  is  customary, 
the  subscript  D  indicates  a  dose  weighted  average.  The  dose  mean  specific  energy  per  event  is 
further  described  and  interpreted  in  section  D1  in  terms  of  probability  density  functions.  It 
should  also  be  noted  from  equation  (6)  that  TD  contains  the  information  describing  the  variability 
of  the  interaction  of  a  single  particle  with  the  site. 

The  purpose  of  this  section  is  to  derive  an  expression  for  the  dose  mean  specific  energy 
per  event  in  terms  of  the  random  physical  variables  described  in  Bl.  The  final  expression  will 
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depend  on  variables  associated  with  the  LET,  path  length  and  collision  magnitude  of  ions  in  the 
site,  which  are  denoted  by  L,  s  and  &,  respectively.  The  distributions  of  these  random  variables 
are  assumed  to  be  in  the  form  of  the  ratio  of  the  second  to  the  first  moment  of  the  respective 
distribution.  Additionally,  it  will  be  shown  how  to  convert  the  description  of  the  fluctuations 
in  the  total  energy  deposited  in  the  site  to  one  which  describes  the  fluctuations  of  the  ionization 
process. 

First  of  all,  consider  the  variability  associated  with  the  LET,  L,  and  path  length,  s, 
variables  for  a  single  ion  event  as  shown  in  figure  lb.  L  and  s  are  independent,  so  the  average 
of  their  product  is  equal  to  the  product  of  their  averages.  Furthermore,  the  product  of  the  two 
variables  is  the  energy  deposited  in  the  site  by  one  ion,  e. 


e  =  L’S  =  L*s 


(38) 


For  similar  reasons, 


s2  =  L2s 2  =  L2s2 


(39) 


Making  use  of  equation  (13),  the  variance  of  e  due  to  L  and  s  is 


Vflr*.,(e)  =  L1  's2  -  (L)2-(s)2 


(40) 


Equation  (40)  gives  the  variance  of  the  energy  deposited  due  to  the  LET  distribution  of  the 
incident  ion  and  its  possible  path  lengths  through  the  volume.  However,  the  energy  deposited 
by  a  single  ion  would  fluctuate  even  if  L  and  s  were  constant.  The  next  step,  therefore,  is  to 
determine  the  variance  contributed  by  these  fluctuations,  which  are  due  to  energy-loss  straggling 
as  discussed  in  section  Bl. 

The  energy  lost  in  single  collisions  by  the  incident  ion  is  described  by  the  random  variable 
6  and  the  number  of  such  collisions  in  the  site  by  the  random  variable  n.  The  process  is  shown 
in  figure  lc.  The  total  energy  deposited  by  the  traversal  of  the  ion  is  therefore  given  by  the  sum 

e  =  6j  +  fi2  +  - •  +  bn  (41) 


where  b,  denotes  the  energy  lost  in  the  i-  of  n  total  collisions.  The  collisions  are  assumed  to 
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occur  independently  so  that  the  6,  are  mutually  independent,  identically  distributed  random 
variables.  This  situation  is  analogous  to  that  expressed  by  equation  (12).  The  calculation  of  the 
average  of  e  and  the  variance  of  e  due  to  energy-loss  straggling,  varstr(t),  in  terms  of  the 
random  variables  6  and  n  is  formally  the  same  as  the  calculations  in  section  B2a  and  will  not  be 
repeated  here.  The  results  are 


S  =  6  ‘ft 


(42) 


and 


varjt)  =  var(b)n  +  (6)2-var(n) 


(43) 


Equations  (42)  and  (43)  are  analogous  to  equations  (18)  and  (27),  respectively.  The  subscript 
”str”  in  (43)  indicates  that  the  variance  results  from  the  straggling  process.  Continuing  the 
analogy,  note  that  n  is  governed  by  Poisson  statistics  so  (43)  can  be  simplified  to  a  result  which 
is  mathematically  similar  to  equation  (30). 


var3 =  n -52 


(44) 


The  notable  feature  of  the  two  processes  described  by  equations  (30)  and  (44)  is  that  they  are 
both  compound  Poisson.  This  results  in  a  similar  mathematical  description  for  the  total  energy 
deposited  in  the  site  and  the  energy-loss  straggling  of  a  single  particle. 

The  variance  expressed  by  (44)  due  to  energy-loss  straggling  occurs  regardless  of  the  LET 
and  path  length  of  the  incident  ion  so  the  processes  are  independent.  Although  this  is  not  a 
sufficient  condition  to  simply  sum  the  (weighted)  variances  of  the  processes,  a  simple  summation 
can  be  justified  for  this  situation.  Further  discussion  of  this  point  is  provided  at  the  beginning 
of  section  Cl.  The  total  variance  of  these  two  processes  is  therefore 


var(e)  =  LZ’S2  -  {L)2is)2  +  n-62 


(45) 


Equation  (45)  represents  the  total  variance  of  the  energy  deposited  in  the  site  by  a  single  particle 
and  includes  the  contributions  due  to  both  ionizations  and  inelastic  excitations. 

If  only  the  ionization  process  is  of  interest,  a  variance  must  be  added  to  equation  (45) 


18 


which  accounts  for  fluctuations  in  the  partitioning  of  the  energy  between  ionizations  and  inelastic 
excitations.  Three  possible  outcomes  of  the  interaction  of  an  incident  ion  with  a  bound  electron 
are  shown  in  figure  Id.  First  of  all,  the  interaction  may  produce  only  excitation.  Secondly,  it 
may  produce  only  a  single  ionization.  Thirdly,  it  may  produce  an  ionization  in  which  the 
resulting  electron  has  kinetic  energy.  The  ejected  electron  subsequently  undergoes  interactions 
which  can  be  ionizations  and/or  excitations.  These  three  outcomes  correspond  to  the  situation 
in  which  the  energy  transferred  to  the  bound  electron  is  less  than,  equal  to  and  greater  than  its 
binding  energy,  respectively.  The  result  of  this  energy  partition  between  ionizations  and 
excitations  is  that  even  for  a  fixed  amount  of  energy  deposited  in  the  site  by  a  single  ion,  the 
number  of  electron-hole  pairs  produced  fluctuate  around  an  average  value.  These  fluctuations 
are  referred  to  as  Fano  fluctuations. 

The  energy  deposited  in  the  site  by  a  single  ion  can  be  written  in  terms  of  the  number  of 
electron-hole  pairs  produced,  q. 


e  «  w 


(46) 


where  w  is  the  average  energy  required  to  produce  one  electron-hole  pair.  In  this  expression, 
c  and  r\  are  treated  as  random  variables  and  w  as  a  constant.  The  average  of  e  is  then  given  by 


e  *  vt) 


(47) 


The  corresponding  variance  of  e  due  to  Fano  fluctuations,  varF(e),  is  then  given  by 


var/e)  =  (w)2-var(Ti) 


(48) 


The  Fano  factor,  F,  is  often  introduced  as  an  empirical  constant  by  which  the  Poisson 
variance  of  x\  must  be  multiplied  to  give  the  observed  variance  of  q.7 


varOl)  =  F-r\ 


OzFz  1 


(49) 


For  F  =  1,  the  process  is  Poisson.  For  F  =  0,  there  is  no  variance,  so  the  same  number  of 
electron-hole  pairs  is  always  produced.  This  implies  that  all  the  energy  is  used  for  ionization. 
Using  equation  (49),  (48)  can  be  rewritten 
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(50) 


var/e)  =  F-q<w)2 


The  Fano  fluctuations  are  assumed  to  be  independent  of  the  fluctuations  due  to  LET,  path  length 
and  energy-loss  straggling  and  the  total  variance  of  the  ionization  process  is  taken  as  the  sum 
of  the  variances  expressed  by  equations  (45)  and  (50). 


vw(e)  =  L2‘S2  -  (L)2is)2  +  tt‘62  +  F-r\{w)2 


(51) 


An  explanation  of  why  these  three  variances  are  simply  summed  is  given  at  the  beginning  of 
section  Cl.  Dividing  both  sides  of  (51)  by  the  average  energy  deposited  by  a  single  ion  and 
expressing  var(t)  in  terms  of  its  first  and  second  moments,  it  is  seen  that 

e2  -  (e)2  =  L2  s2  _  (L)2js)2  +  n&2  +  F-i\jw)2  (52) 

e  e  e  e  e 

Note  from  (38)  that  the  second  terms  on  each  side  of  the  above  equation  are  identical. 
Therefore, 


e2  _  L2 • s 2  +  n-62  +  F •  r| -(w)2 
e  e  e  e 


(53) 


Substituting  the  expressions  given  by  (38),  (42)  and  (47)  for  the  average  of  e  into  the  first, 
second  and  third  terms,  respectively,  of  the  right  hand  side  of  (53)  results  in 


E 


L  s 


5 


+  F'W 


(54) 


Equation  (54)  expresses  the  moment  ratio  of  the  energy  deposited  in  the  site  by  single  ions  in 
terms  of  the  known  moment  ratios  of  the  LET,  path  length  and  collision  magnitude  distributions, 
and  known  parameters  describing  the  Fano  fluctuations.  As  discussed  previously  in  connection 
with  equation  (10),  the  various  moment  ratios  in  equation  (54)  can  be  interpreted  as  weighted 
averages  of  the  appropriate  distribution.  Combining  this  result  with  equation  (36)  gives  the  final 
expression  for  the  dose  mean  specific  energy  per  event. 
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—  1.6xlO'1J 

=  KfPo 


‘  1}  j2  »2  _  - 

— —  +  —  +  F-w 


V  L  s 


(55) 


Equations  (35)  and  (55)  give  the  relative  variance  of  the  dose  deposited  in  the  site  in 
terms  of  known  or  calculable  quantities.  Equation  (55)  characterizes  the  variability  of  a  single 
ion  event  in  the  site  while  the  dependence  of  the  relative  variance  on  the  average  dose  is 
contained  entirely  in  (35).  These  two  equations  describe  the  fluctuations  in  the  ionization  process 
(electron-hole  pair  production)  which  occur  in  the  small  volume.  To  describe  the  relative 
variance  due  to  the  complete  energy  deposition  process,  which  would  include  both  ionizations 
and  inelastic  excitations,  the  same  two  equations  are  used  except  that  the  term  F-w  is  set  equal 
to  zero  in  (55).  The  results  derived  above  are  equivalent  to  those  presented  by  Kellerer.1,2 

Finally,  a  note  about  units  in  equation  (55).  The  ratio  of  the  second  moment  to  the  first 
moment  of  the  incident  particle’s  LET  is  in  units  of  keV/pm,  the  ratio  of  the  path  length 
moments  through  the  site  is  in  pm,  the  ratio  of  the  moments  of  the  distribution  of  collision 
magnitudes  is  in  keV,  F  is  dimensionless,  w  is  in  keV,  V0  is  in  cm3  and  p0  is  in  g/cm3.  The 
factor  of  1.6xl0*11  simply  converts  the  absorbed  dose  from  keV/g  to  rads.  The  dose  mean 
specific  energy  per  event  is  in  rad  units. 


3.  Energy  Deposition  Fluctuations  with  Random  Variables  Described  by  Relative  Variances 

For  some  applications  it  is  more  convenient  to  describe  the  known  distributions  of  the 
random  variables  in  terms  of  relative  variances  rather  than  in  terms  of  moment  ratios.  One  of 
the  advantages  of  using  relative  variances  is  that  they  are  dimensionless.  In  section  3a,  relations 
between  the  moment  ratios  and  relative  variances  of  the  pertinent  random  variables  are  obtained. 
Then,  using  these  relations  in  section  3b,  equations  (35)  and  (55)  are  converted  into  a  form  in 
which  the  variability  of  the  dose  deposited  in  the  site  is  expressed  in  terms  of  relative  variances. 


a.  Relations  between  Moment  Ratios  and  Relative  Variances  of  the  Random 
Variables 

In  this  section,  relations  between  the  ratio  of  the  second  moment  to  the  first  moment  and 
the  relative  variance  are  obtained  for  the  LET,  path  length  and  energy-loss  straggling 
distributions.  Additionally,  the  term  corresponding  to  the  Fano  fluctuations  in  equation  (55)  is 
also  expressed  in  terms  of  its  relative  variance. 

Equation  (6)  gives  the  general  relation  between  the  relative  variance  and  the  moment  ratio 
for  a  continuous  random  variable.  Applying  it  to  the  variables  L  and  s  results  in 


21 


(56) 


~  =  L<yL+ 1) 

L 

=  s  +  1)  (57) 


Note  that  when  the  variance  of  L  or  s  is  zero,  the  corresponding  moment  ratio  is  simply  the 
average  value.  These  limiting  cases  correspond  to  the  situation  of  a  monoenergetic  source  of 
ions,  and  to  the  situation  where  all  ion  paths  through  the  site  have  the  same  length. 

As  discussed  in  Bl,  the  energy-loss  straggling  process  itself  is  a  compound  Poisson 
process  because  the  ion  undergoes  n  independent  collisions  of  varying  magnitude  6  at  each 
collision.  The  variance  of  the  energy-loss  straggling  process,  given  by  (44),  was  obtained  by 
using  Burgess’  variance  theorem  to  combine  the  random  variables  n  and  6.  The  relative 
variance,  Vttr,  is  obtained  by  dividing  the  variance  by  the  square  of  the  average.  Since  the 
average  is  given  by  (42),  the  relative  variance  due  to  energy-loss  straggling  is 


(58) 


The  explicit  dependence  on  the  average  number  of  collisions  can  be  eliminated  using  (42).  The 
relation  between  the  moment  ratio  of  the  distribution  of  collision  magnitudes  and  the  relative 
variance  of  the  straggling  process  is  therefore 


& 

5 


(59) 


Although  there  is  no  explicit  dependence  on  n  in  equation  (59),  it  is  understood  that  Vjrr  is  the 
relative  variance  for  the  entire  energy-loss  straggling  process,  which  includes  contributions  due 
to  fluctuations  in  both  the  number  and  magnitude  of  collisions  as  the  ion  crosses  the  site. 
Although  equation  (59)  is  a  more  useful  relation,  it  is  interesting  to  write  Vstr  in  terms  of  the 
relative  variance  of  the  collision  magnitudes,  V6.  Applying  equation  (6)  to  the  variable  6  and 
substituting  the  result  into  (58),  the  following  equation  is  obtained. 
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v„  -  i<i  *  vt) 

n 


(60) 


Since  n  is  governed  by  Poisson  statistics,  the  first  term  in  (60)  is  the  relative  variance  of  the 
number  of  collisions.  The  second  term  therefore  describes  the  variability  due  to  the  collision 
magnitudes. 

Finally,  the  term  involving  the  Fano  fluctuations  in  equation  (55)  will  now  be  related  to 
the  relative  variance  of  the  Fano  fluctuations,  VF.  Equation  (50)  gives  the  variance  of  the 
fluctuations  while  (47)  gives  the  average.  Again,  the  relative  variance  is  obtained  by  dividing 
the  variance  by  the  square  of  the  average. 


VF  =  4  (61) 

11 

This  would  also  follow  directly  from  equation  (49).  Using  (47)  to  eliminate  the  dependence  on 
the  average  number  of  electron-hole  pairs,  rf,  the  following  result  is  obtained. 

F 'W  -  ~t'VF  (62) 


b.  Conversion  to  the  Relative  Variance  Description 

The  relations  necessary  to  convert  the  description  with  the  random  variables  expressed 
as  moment  ratios  (equations  35  and  55)  to  one  with  the  random  variables  expressed  as  relative 
variances  are  given  by  (56),  (57),  (59)  and  (62).  Substituting  these  relations  into  equation  (55) 
and  with  some  algebraic  manipulation  it  is  seen  that 


Zn  -  1-6.X-1Q;11  \L~siVL  +  l)iVs  +  1)  ♦  e  -Vm  +  fVF] 


VPo 


(63) 


Noting  from  (38)  that  the  average  LET  times  the  average  path  length  is  the  average  energy 
deposited  by  a  single  ion,  T,  this  can  be  written  as 
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•«-(i +  vL  +  k  +  vlK *  ym  +  yF) 


(64) 


= 


1.6xlQ-n 

VPo 


Recall  from  (33)  that  the  multiplicative  factor  l.6x  10'11/Vo-p0  converts  energy  deposition  in  keV 
to  dose  in  rads.  Equation  (64)  can  then  be  expressed  as 


zd  =  a*(i  +  vi  +  v.  +  v*K +  y. +  yF) 


(65) 


where  d  is  the  average  dose  deposited  in  the  site  by  a  single  event. 

The  sum  of  all  the  relative  variance  terms  inside  the  parenthesis  on  the  right  hand  side 
of  (65)  is  called  the  single  event  relative  variance,  V{. 

Vl-Vl*V.*  VL  V,  *V„*VF  <«® 


This  quantity  is  the  relative  variance  of  the  ionization  produced  by  a  single  particle  in  the  small 
volume,  excluding  that  associated  with  Poisson  statistics.  It  is  therefore  independent  of  dose. 
Writing  equation  (65)  in  terms  of  V,  gives  the  useful  result 


zD  =  dil  +  Vl) 


(67) 


The  next  step  is  to  express  the  relative  variance  of  the  dose  deposited  in  the  site  in  terms 
of  Vj.  To  accomplish  this,  (67)  is  substituted  into  (35)  to  obtain 


(68) 


Since  d  is  the  dose  deposited  by  a  single  particle  and  D  is  the  total  dose,  this  can  be  rewritten 
as 
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VD  =  i-d  +  iy 


(69) 


where  N  is  the  average  number  of  incident  ions  which  strike  the  site.  The  mathematical 
similarity  of  equations  (60)  and  (69)  is  apparent.  Again,  this  results  from  the  fact  that  the 
energy-loss  straggling  of  an  ion  crossing  a  site  and  the  deposition  of  dose  in  a  site  are  both 
compound  Poisson  processes. 

Equations  (66)  and  (69)  are  the  final  results.  They  relate  the  relative  variance  of  the  dose 
deposited  in  the  site  to  known  relative  variances  and  the  average  number  of  incident  particles 
which  strike  the  site.  The  latter  quantity  can  be  calculated  from  the  average  deposited  dose. 
Equation  (66)  characterizes  the  variability  of  a  single  particle  interaction  in  the  site  and  (69) 
relates  the  total  relative  variance  to  the  single  event  relative  variance  and  the  number  of  events. 
All  quantities  in  these  two  equations  are  dimensionless.  Again,  these  results  are  equivalent  to 
those  presented  by  Kellerer.1,2 

Equations  (66)  and  (69),  with  the  Fano  fluctuations  accounted  for  by  VF,  describe  the 
fluctuations  in  the  ionization  process  (electron-hole  pair  production).  To  describe  the  fluctuations 
due  to  both  ionizations  and  inelastic  excitations,  the  same  two  equations  are  used  except  that  VF 
is  set  equal  to  zero  in  (66). 

It  is  interesting  to  note  two  other  features  of  this  description  using  relative  variances. 
First,  setting  Vj  =  0  in  (69),  which  is  equivalent  to  setting  the  relative  variances  due  to  LET, 
path  length,  straggling  and  Fano  fluctuations  equal  to  zero,  implies  that  the  relative  variance  of 
the  dose  deposited  in  the  site  is  described  entirely  by  Poisson  statistics.  Therefore,  the  first  term 
in  (69)  is  identified  with  Poisson  fluctuations  in  the  number  of  particles  that  strike  the  site. 

Secondly,  note  that  the  variability  due  to  L  and  s,  which  occurs  as  a  result  of  the  product 
of  the  two  independent  random  variables,  results  in  three  terms  involving  their  relative  variances 
in  equation  (66).  This  is  the  general  result  for  the  product  of  two  independent  random 
variables.1,8  It  is  usually  assumed  in  texts  on  statistics  that  the  relative  variances  are  small  and 
an  expression  is  derived  from  a  series  expansion  which  neglects  the  term  involving  the  product 
of  the  two  relative  variances.  The  expression  presented  here  is  more  general  and  includes  the 
possibility  that  either  or  both  of  VL  and  Vs  can  be  large. 

Finally,  note  that  equation  (67)  is  a  convenient  expression  for  converting  between  the 
relative  variance  description  given  by  equations  (66)  and  (69)  and  the  moment  ratio  description 
given  by  equations  (35)  and  (55).  It  allows  the  dose  mean  specific  energy  per  event  to  be 
calculated  given  the  single  event  relative  variance,  V,,  or  its  components  VL,  v„  Vstr  and  Vf. 
Conversely,  it  allows  Vx  to  be  calculated  given  the  dose  mean  specific  energy  per  event  or  the 
moment  ratios  of  the  random  variables  L,  s  and  6,  and  the  product  of  the  Fano  factor  and  the 
average  energy  to  produce  an  electron-hole  pair. 
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C.  The  Two  Component  Model 


The  relations  in  section  B  describe  energy  deposition  fluctuations  for  a  single  type  of 
event  such  as  the  traversal  of  a  site  by  ions.  These  relations  were  published  at  a  time  when 
interest  in  radiobiology  was  centered  on  the  micrometer  dimensions  of  cell  nuclei.  For  sub- 
micrometer  dimensions,  however,  it  is  no  longer  sufficient  to  consider  only  the  incident  ions 
which  directly  strike  the  site.9,10  The  range  of  secondary  electrons  in  the  medium  can  be 
significant  compared  to  the  site  dimensions,  and  energy  transfer  by  electrons  into  and  out  of  the 
small  volume  must  be  accounted  for.  In  some  cases  this  contribution  becomes  equally  important 
to  that  of  the  direct  ion  interactions  in  the  site.  An  electron  event  (also  commonly  referred  to 
as  an  indirect  event  or  a  "toucher")  is  considered  to  occur  when  a  secondary  electron  strikes  the 
target  site  after  being  created  by  an  incident  ion  outside  of  the  site.  By  consideration  of  both 
electron  and  ion  events,  an  approach  will  be  developed  to  describe  energy  deposition  fluctuations 
induced  by  ions  in  target  sites  having  dimensions  as  small  as  one  nanometer. 

Section  1  gives  a  review  of  the  appropriate  probability  theory  which  will  be  applied  to 
combine  the  dose  distribution  of  ions  that  directly  interact  with  the  site  and  the  dose  distribution 
of  secondary  electrons  that  interact  with  the  site  after  being  produced  by  ions  in  the  surrounding 
medium.  The  combination  of  the  two  distributions  is  a  key  part  of  the  derivation  which  has  not 
been  previously  addressed  in  the  context  of  microdosimetry.  In  section  2,  the  Two  Component 
Model  is  obtained  for  the  situation  where  the  random  variables  of  the  energy  deposition  process 
are  described  by  moments  of  their  respective  distributions.  Section  3  describes  the  Two 
Component  Model  where  the  random  variables  are  described  by  relative  variances. 


1.  Combined  or  Mixed  Probability  Distributions 

The  Two  Component  Model  is  based,  from  a  statistical  viewpoint,  on  combining 
parameters  such  as  the  mean  and  variance  of  two  independent  probability  distributions  to  obtain 
the  corresponding  parameters  of  the  combined  or  mixed  distribution.  It  will  be  necessary  to 
distinguish  between  the  situation  where  the  following  analysis  of  combined  probability 
distributions  is  applicable  and  when  it  is  possible  to  simply  sum  the  weighted  variances  of  the 
component  distributions.  Breitenberger11  makes  a  clear  distinction  between  the  two  cases.  The 
analysis  of  this  section  applies  to  what  he  calls  "alternative  events".  For  alternative  events,  one 
and  only  one  of  two  or  more  possible  alternatives  will  occur  in  a  given  trial.  As  an  example  in 
the  context  of  this  work,  consider  the  dose  deposited  in  a  target  site  resulting  from  ion  irradiation 
of  the  material.  An  ion  may  directly  strike  the  site,  and  this  process  has  an  associated  dose 
distribution  per  event.  Alternatively,  the  ion  may  miss  the  site  but  create  a  secondary  electron 
which  deposits  dose  in  the  site.  The  two  independent  dose  distributions  due  to  ions  and  electrons 
result  from  alternative  events,  and  the  variance  of  the  combined  distribution  is  not  simply  the 
sum  of  the  weighted  variances  of  the  component  distributions. 

On  the  other  hand,  it  is  possible  to  simply  sum  the  weighted  variances  for  the  situation 
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A -distribution: 
Na  events 
mean  xA 
variance  var(xA) 


B- distribution: 
Nb  events 
mean  xB 
variance  vartxg) 
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that  Breitenberger  calls  "simultaneous  events".  Simultaneous  events  are  independent  events 
associated  with  the  same  ion  that  happen  at  the  same  time.  For  example,  consider  the 
fluctuations  in  electron  -  hole  pair  production  resulting  from  a  single  ion  event  as  discussed  in 
section  B2b.  The  three  variances  which  contribute  to  the  fluctuations  are  due  to  the  LET  times 
the  path  length  of  the  ion  [equation  (40)],  energy-loss  straggling  [equation  (44)]  and  Fano 
fluctuations  [equation  (50)].  These  three  processes  occur  simultaneously  and  independently  as 
the  ion  crosses  the  site.  Under  these  conditions,  the  total  variance  is  given  by  (51),  which  is  the 
sum  of  the  three  individual  variances. 

The  following  analysis  is  for  a  combined  probability  distribution  of  alternative  events. 
Consider  two  independent  distributions  of  continuous  random  variables  xA  and  xB  shown  in  figure 
2.  For  clarity,  they  are  shown  separated,  but  they  may  be  overlapping.  They  will  be  referred 
to  as  the  /I -distribution  and  the  B-distribution.  The  -distribution  is  made  up  of  NA  events  which 
have  a  mean  value  of  x\  and  a  variance  of  var(xj.  Similarly,  the  B-distribution  consists  of  NB 
events  having  a  mean  xj  and  a  variance  varft*).  The  ordinate  shown  in  the  figure  shows  the 
probability  density  function, /ft).  It  represents  the  number  of  events  which  have  a  value  of  x, 
normalized  to  the  total  number  of  events,  NA  +  NB.  In  other  words,  it  is  the  fraction  of  the  total 
number  of  events  which  have  a  value  of  x.  This  type  of  probability  function  is  often  called  a 
frequency  weighted  probability  density  function  in  the  context  of  microdosimetry  theory  and  will 
be  elaborated  on  in  section  D.  It  is  assumed  the  means,  variances  and  number  of  events 
associated  with  each  component  distribution  are  known.  The  goal  is  to  combine  these  known 
parameters  of  the  two  component  distributions  to  arrive  at  the  mean  and  variance  of  the  overall 
mixed  distribution.  Equivalently,  the  goal  is  to  combine  the  moments  of  the  component 
distributions  to  obtain  the  moments  of  the  mixed  distribution. 

The  problem  of  combined  or  mixed  distributions  has  been  treated  in  the  context  of  applied 
probability  theory  by  Deming.12  The  basic  result  for  combining  moments  of  the  two  component 
distributions  is 


*"•  =  ?•*;  +  (i  -pyx; 


(70) 


where  xAm  and  xBm  are  the  m-  moments  of  the  A-  and  B-distributions,  respectively,  and  x"  is  the 
m-  moment  of  the  mixed  distribution.  The  weighting  factors  P  and  1  -  P  are  the  fraction  of  the 
total  number  of  events  in  the  A-  and  B-distributions,  respectively. 


P  = 


Na  +  Nb 


(71) 
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(72) 


t 


1  -  P  * 


N, 


B 


Na+Nb 


Deming  refers  to  P  and  1  -  P  as  the  normalized  or  relative  frequencies  of  the  A-  and  B- 
distributions. 

The  average,  x,  and  the  variance,  var(x),  of  the  combined  distribution  can  be  obtained 
making  use  of  equation  (70).  The  average  results  from  setting  m  equal  to  1. 


x  *  P  -XA  +  (X-P)  XB  (73) 


where  ~A  and  x~B  are  the  averages  of  the  A-  and  B-distributions,  respectively.  The  variance  of  the 
combined  distribution  can  be  obtained  by  starting  with  equation  (4)  and  inserting  expressions  for 
the  first  and  second  moments  of  x  obtained  using  equation  (70).  This  gives 

var{x)  =  P  -xA  +  (1  -  P)  x*  -  [P  •  TA  +  (1  -  P)-Tg  f  <74) 

The  first  two  terms  on  the  right-hand  side  of  equation  (74)  can  be  expressed  in  terms  of  the 
average  and  variance  of  the  appropriate  distribution  by  making  use  of  (4). 


var(x)  =  P\var{xA)  -  (xf\  +  (1  -  P)[var(xJ  -  (xj2]  -  [P  xA  *  (1  -  P)  J^  (75) 


Note  that  the  variance  of  the  mixed  distribution  is  now  written  in  terms  of  the  averages  and 
variances  of  the  component  distributions.  All  that  remains  is  to  simplify  equation  (75). 
Expanding  the  last  term  and  performing  some  algebraic  manipulations,  it  can  be  shown  that 

var(x)  =  P  •  var(xA)  +  (1  -  P)-\arix^  +  P  (1  -  P)-(xg  -  x^)2  (76) 


This  result,  along  with  alternative  expressions,  is  given  by  Deming.  In  addition,  expressions  are 
given  for  the  variance  of  an  arbitrary  number  of  component  distributions.  In  equation  (76), 
varfxj  and  var(Xg)  are  referred  to  as  the  "internal"  variances  and  the  last  term  is  called  the 
"external"  variance.  The  variance  of  the  combined  distribution  is  the  sum  of  the  internal 
variances,  properly  weighted,  and  the  external  variance,  which  is  proportional  to  the  square  of 
the  separation  of  the  component  averages. 
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In  summary,  for  a  combined  or  mixed  distribution  of  alternative  events,  equations  (70), 
(73)  and  (76)  give  the  moments,  mean  and  variance  of  the  combined  distribution  in  terms  of 
known  quantities  of  the  independent  component  distributions.  The  weighting  factors  P  and  1  - 
P  are  given  by  (71)  and  (72),  respectively.  It  has  been  shown  that  the  correct  procedure  for  this 
situation  is  to  sum  the  weighted  moments  of  the  component  distributions,  not  the  weighted 
variances.  Simply  summing  the  weighted  variances  is  a  reasonable  approximation  only  if  the 
averages  of  the  component  distributions  are  nearly  equal. 


2.  Energy  Deposition  Fluctuations  with  Random  Variables  Described  by  Moments 

The  new  microdosimetry  model  explicitly  accounts  for  two  components  which  contribute 
to  the  fluctuations  in  energy  deposited  in  the  site.  The  components  are  shown  schematically  in 
Figure  3.  The  first  contribution  is  due  to  incident  ions  which  directly  strike  the  site.  This  will 
subsequently  be  referred  to  as  the  ion  or  direct  component  and  denoted  with  the  subscript  "ion". 
As  shown  in  the  figure,  some  of  the  energy  initially  deposited  in  the  site  may  be  carried  out  of 
the  site  by  secondary  electrons.  The  second  contribution  is  due  to  secondary  electrons  which 
strike  the  site  after  being  created  by  ions  outside  the  site.  This  will  be  referred  to  as  the  electron 
or  indirect  component  and  indicated  by  the  subscript  “e".  This  nomenclature  and  notation  is 
similar  to  that  used  previously  by  Kellerer,13  the  present  authors14  and  Olko  and  Booz.ls  Both 
the  ion  and  the  electron  components  contribute  to  the  fluctuations  in  energy  deposited  in  the 
target  site.  In  order  to  characterize  the  overall  energy  distribution,  the  means  and  variances  of 
the  ion  and  electron  distributions  must  be  properly  combined.  The  mathematical  problem  of 
combining  the  two  probability  distributions  has  been  discussed  in  the  previous  section. 

To  begin,  note  that  (69)  expresses  the  total  relative  variance  of  the  dose  in  terms  of  the 
relative  variance  due  to  single  events  and  the  average  number  of  events  in  the  site.  Using  the 
definition  of  relative  variance  given  by  (5),  equation  (69)  can  be  written 


var(D )  _  _1_ 

(D)2  N 

where  D  and  var(D)  are  the  mean  and  variance  of  the  total  dose  deposited  in  the  site,  and  <Tand 
var(d)  are  the  mean  and  variance  of  the  dose  deposited  by  a  single  particle. 

Furthermore,  equation  (18)  can  be  expressed  in  terms  of  dose  related  quantities. 

D  =  d-N  (78) 


vaijd) 


W 


(77) 


Inserting  (78)  into  (77),  the  variance  of  the  dose  can  be  written  as 
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Figure  3.  Factors  contributing  to  energy  deposition  in  the  Two  Component  Model. 


varfD)  =  N  •  [(d)1  +  var(d)\ 


(79) 


Since  the  variance  of  d  is  given  by 


var(d)  -  d1  -  (d)2 


(80) 


equation  (79)  simplifies  to 


var(D)  *  N-d2 


(81) 


The  variance  of  the  total  dose  deposited  in  the  site  is  the  product  of  the  average  number  of  events 
and  the  second  moment  of  the  single  event  dose  distribution. 

The  problem  is  to  determine  an  expression  for  the  total  variance  for  the  situation  where 
the  two  multiplicative  factors  in  (81)  result  from  a  combined  distribution  of  ion  and  electron 
events.  The  first  factor,  the  average  number  of  total  events,  is  the  sum  of  the  average  number 
of  ion  events  and  the  average  number  of  electron  events. 

N  =  +  We  (82) 


The  remaining  factor,  the  second  moment  of  the  single  event  dose  distribution,  is  now 
considered.  Any  given  event,  by  definition,  deposits  dose  in  the  site.  The  dose  deposited  by 
this  event  results  from  either  an  ion  striking  the  site  qt  an  electron  striking  the  site.  There  are 
two  independent  dose  distributions  which  result  from  the  two  types  of  events.  Therefore,  as 
discussed  in  section  Cl,  the  ion  and  electron  events  are  "alternative"  events.  Equation  (70)  can 
then  be  applied  to  combine  the  moments  of  the  ion  and  electron  single  event  dose 
distributions  to  obtain  the  nP-  moment  of  the  overall  dose  distribution.  The  result  for  combining 
the  second  moments  is 


d1  =  />•<&  *  (1  -P)  d] 


(83) 


where  dion  and  de  are  the  dose  deposited  in  the  site  by  single  ion  and  by  single  electron  events, 
respectively.  P  and  1  -  P  represent  the  fractions  of  the  total  number  of  events  due  to  ions  and 
electrons,  respectively.  Equation  (81)  can  now  be  written  as 

vaiiD )  =  N  •[?  • dl,  +  (1  -  P) 'd]\  (84) 

It  turns  out  to  be  more  convenient  to  express  the  results  in  terms  of  dose  fractions  rather 
than  the  event  fractions  P  and  1  -  P.  Next,  expressions  will  be  obtained  which  relate  the  two 
types  of  fractions.  The  fraction  of  direct  ion  hits  of  the  site  is,  on  the  average,  given  by 


P  = 


(85) 


The  corresponding  fraction  of  electron  hits  is 


1 


-  P  = 


+ 


(86) 


The  dose  fractions  are  given  by  similar  expressions.  If  Dion  is  the  dose  deposited  by  ion  strikes 
and  De  is  the  dose  deposited  by  e’ectron  strikes, 


y*  _  _ ton 

ion  *"  ■■  *=r- 

D.  +  D 

ton  t 


(87) 
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(88) 


D,  +  D 

ion  t 


where^  and  1  -/^  are  the  fractions  of  dose  deposited  by  the  ion  and  electron  components, 
respectively.  Furthermore,  if  dion  is  the  dose  deposited  per  ion  event  and  de  is  the  dose  deposited 
per  electron  event,  then 


Djon  =  dion  ion 


Dt  =  dt  • Nt 


(89) 


(90) 


Substituting  (89)  and  (90)  into  (87)  and  also  into  (88)  results  in 


/*.- 


dion  ’  Nion 


d,  *  N,  +  d  ‘N 

ion  ion  e  e 


(91) 


dm -N. 

1  ~,km  *  +  K 


(92) 


Dividing  the  numerator  and  denominator  of  these  two  equations  by  Nion  +  Ne,  they  become 


,  m  _ Pdion 

J  ion  _  ,  ■ 


P -dLMl-fK 


(93) 


(1-P)<*, 

?  a;  +  a -F)Jt 


(94) 
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The  denominator  of  the  above  two  expressions  is  the  average  value  of  the  dose  deposited  per 
event  of  the  combined  distribution.  This  can  be  seen  from  equation  (73),  the  general  expression 
for  the  average  of  the  mixed  distribution. 


d  =  P -<k  +  (l-fK 


(95) 


Using  (95),  equations  (93)  and  (94)  can  be  written 


fion  ~ 


Pd 


ion 


(96) 


1  -fion* 


(1  -P)de 


(97) 


These  last  two  expressions  are  the  relations  between  the  event  fractions,  P  and  1  -  P,  and  the 
dose  fractions,  fm  and  1  -  fion.  Solving  (96)  for  P  and  (97)  for  1  -  P,  and  inserting  the  results 
into  (84),  the  variance  of  D  becomes 


var(D)  =  N- 


+ 


a  -/j' d 
dt 


(98) 


Equation  (78)  can  be  used  to  write  this  in  terms  of  the  average  dose,  which  is  a  measurable 
quantity. 


-dr  dr 


(99) 


*ion 


Recall  from  (37)  that  the  dose  mean  specific  energy  per  event  is  the  ratio  of  the  second  moment 
to  the  first  moment  of  the  dose  deposited  by  a  single  particle  in  the  site.  It  is  natural,  then,  to 
identify  the  two  moment  ratios  in  the  above  equation  as  analogous  quantities  for  ions  and 
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electrons. 


(100) 


-  < 
ZD*  ~  V 


(101) 


In  microdosimetry  terminology,  equation  (100)  represents  the  dose  mean  specific  energy  per 
event  for  direct  ion  strikes  of  the  site,  and  (101)  is  the  analogous  quantity  for  electron  strikes. 
These  two  quantities  are  further  described  and  interpreted  in  section  D. 

The  variance  given  by  (99)  can  now  be  written 


Var<D)  =fu»‘D'ZDjo*  +  (!  ~fiJ  D  ZD* 


(102) 


This  can  be  converted  to  a  relative  variance  by  dividing  both  sides  of  (102)  by  the  square  of  the 
average  deposited  dose. 


fjon  ZDjon  *  (1  fjot^  ZDjc 

D 


(103) 


This  result  can  be  compared  to  equation  (35).  It  is  seen  that  in  order  to  obtain  the  total  relative 
variance  of  the  dose,  the  dose  mean  specific  energies  of  ions  and  electrons  must  be  weighted  by 
the  fraction  of  dose  deposited  in  the  site  by  each  type  of  event. 

When  fon  is  close  to  unity,  the  ion  component  dominates.  For  this  situation  nearly  all  the 
energy  initially  deposited  in  the  site  remains  within  it.  As  the  energy  transfer  by  secondary 
electrons  into  and  out  of  the  site  becomes  significant, becomes  less  than  one,  and  the  electron 
component  becomes  more  important.  Intuitively,  it  would  be  expected  that  the  electron 
component  would  be  most  important  for  small  site  dimensions  and  for  high  energy  incident  ions. 

The  dose  mean  specific  energies  per  event  for  ions  and  electrons  are  obtained  by  applying 
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equation  (55)  to  each  type  of  event. 
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(104) 


(105) 


where  Lion  and  are  the  random  variables  which  represent  the  LET  of  the  incident  ion  and 
incident  secondary  electron,  respectively,  sion  and  se  are  the  path  lengths  of  the  ion  and  electron 
in  the  site,  6ion  and  6e  represent  the  energy  deposited  in  a  single  ion  collision  and  a  single 
electron  collision  in  the  site,  Fim  and  Fe  are  the  Fano  factors  for  the  two  components,  w~  and 
We  are  the  average  energies  required  to  produce  one  electron-hole  pair  by  each  component,  V0 
is  the  volume  of  the  site  and  p0  its  density.  The  presence  of  the  factor  flon  in  equation  (104) 
accounts  for  transport  of  energy  initially  within  the  site  to  the  outside,  as  indicated  in  figure  3. 
An  analogous  factor  is  not  present  in  equation  (105)  because  the  overwhelming  majority  of 
incident  secondary  electrons  do  not  create  tertiary  electrons  within  the  site  with  sufficient  energy 
to  escape.16 

The  last  term  in  each  of  (104)  and  (105)  is  associated  with  the  Fano  fluctuations  of  the 
ion  and  electron  components,  respectively.  Setting  them  equal  to  zero  results  in  a  description 
of  the  energy  deposition  (ionizations  and  excitations)  in  the  site  while  including  them  describes 
the  ionization  (electron-hole  pair  production)  in  the  site.  The  Fano  fluctuations  are  a  very  small 
contribution  for  micrometer-sized  sites.2  However,  they  can  become  significant  as  site 
dimensions  approach  nanometers. 

Equations  (103),  (104)  and  (105)  are  the  basic  equations  for  the  Two  Component  Model 
when  the  distributions  associated  with  the  random  variables  for  the  ion  and  electron  components 
are  given  in  terms  of  moment  ratios.  The  latter  two  equations  describe  the  interaction  of  single 
particles  (ions  and  electrons,  respectively)  with  the  site.  Equation  (103)  weights  the  two 
components  properly  and  accounts  for  the  average  dose  dependence  of  the  fluctuations. 
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3.  Energy  Deposition  Fluctuations  with  Random  Variables  Described  by  Relative  Variances 

In  the  discussion  of  section  B,  it  was  noted  that  the  mathematical  manipulation  may  be 
easier  if  the  contributions  to  the  fluctuations  are  described  in  terms  of  relative  variances  rather 
than  moment  ratios.  In  this  section,  the  moment  ratio  description  of  the  previous  section  is 
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converted  to  a  form  in  which  the  random  variables  associated  with  the  Two  Component  Model 
are  given  as  relative  variances.  Note  that  upon  inspection  of  equation  (103),  the  total  relative 
variance  of  the  dose  deposited  is  the  sum  of  the  relative  variances  due  to  the  ion  and  electron 
components.  In  transforming  to  the  relative  variance  description,  it  will  be  simpler  to  consider 
each  component  relative  variance  separately.  Denoting  them  by  VD  ion  and  VD  e, 


■tr  _  f ion  ^D^on 

v ***  =  —=— 


(106) 
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D 


(107) 


To  convert  the  Two  Component  Model  to  a  description  based  on  relative  variances,  equation  (67) 
is  a  convenient  starting  point.  Applying  it  to  the  ion  component, 


^DJon  ^ ion  +  Joit) 


(108) 


where  dion  is  the  dose  deposited  in  the  site  by  a  single  ion,  allowing  for  transport  of  energy 
initially  within  the  site  to  the  outside,  and  Vi  ion  is  the  corresponding  relative  variance  due  to 
these  single  events.  In  analogy  to  (66),  Vl  ion  is  given  by 

VUm  =  VUon  ♦  Kjon  +  +  VmJm  +  V9Jm  <109> 


where  VL  lon,  Vs  ion,  Vstr  lon  and  VF  ion  are  the  relative  variances  of  LET,  path  length,  straggling  and 
Fano  fluctuations  for  die  incident  ions.  Substituting  (108)  into  (106)  gives  the  following 
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In  analogy  to  equation  (69),  the  factor  in  front  of  the  parenthesis  will  be  expressed  in  terms  of 
the  average  number  of  ion  events.  This  quantity  is  equal  to  the  fraction  of  the  average  dose 
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deposited  by  ions  divided  by  the  average  dose  deposited  per  ion  event. 


(Ill) 


Combining  (110)  and  (111)  results  in 
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(112) 


This  expression  gives  the  relative  variance  of  the  dose  deposited  in  the  site  by  ions  in  terms  of 
dimensionless  quantities. 

A  similar  procedure  can  be  used  to  obtain  the  relative  variance  of  the  dose  deposited  in 
the  site  by  secondary  electrons.  Applying  (67)  to  the  secondary  electrons 


i*»y 


(113) 


where  de  is  the  dose  deposited  in  the  site  by  a  single  electron  created  outside  of  it.  Vl  e  is  the 
relative  variance  due  to  these  single  electron  events  and  is  given  by 

yu  -  yu  +  k,  *  vu  K,  *  ymt  *  y„  <n4» 


VL  e,  Vs  e,  Vstr  e  and  VFe  are  the  relative  variances  of  the  LET,  path  length,  straggling  and  Fano 
fluctuations  of  the  incident  electrons.  Substituting  (113)  into  (107)  results  in 


(115) 
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The  average  number  of  electron  events  is  given  by 


W_(LMR 
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(116) 


Equation  (115)  can  now  be  rewritten  in  terms  of  dimensionless  quantities  by  making  use  of 

(116). 


(117) 


This  gives  the  relative  variance  of  the  dose  deposited  in  the  site  by  the  electron  component.  The 
two  terms  of  equation  (103)  associated  with  the  relative  variance  of  the  ion  and  electron 
components,  properly  weighted,  are  now  given  by  (112)  and  (117),  respectively.  The  total 
relative  variance  is  therefore  given  by  their  sum. 
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Equations  (109),  (114)  and  (118)  are  the  basic  relations  for  evaluating  the  fluctuations  of 
electron-hole  pair  production  in  the  site,  given  that  the  basic  information  of  the  random  variables 
is  in  the  form  of  relative  variances.  Equations  (109)  and  (114)  account  for  the  variability  of  the 
ionizations  in  the  site  produced  by  single  incident  ions  and  single  incident  electrons,  respectively. 
Equation  (118)  properly  weights  the  two  components  and  incorporates  the  average  dose 
dependence  in  terms  of  the  number  of  ion  and  electron  events.  If  energy  fluctuations  are  the 
quantity  of  interest  instead  of  ionizations,  VFion  and  VFe  are  set  equal  to  zero  in  (109)  and  (114), 
respectively. 

Equations  (109),  (114)  and  (118)  are  particularly  well  suited  for  evaluating  the  relative 
contributions  of  the  various  physical  factors  which  contribute  to  the  energy  deposition  (or 
electron-hole  pair)  fluctuations.  The  total  relative  variance  can  be  broken  into  contributions  due 
to  ion  and  electron  events.  Furthermore,  the  contributions  due  to  the  LET  distribution,  path 
length  distribution,  straggling  distribution,  Fano  fluctuations  and  Poisson  statistics  for  both  ion 
and  electron  events  can  be  analyzed. 
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D.  Microdosimetric  Quantities  and  the  Two  Component  Model 


The  microdosimetry  literature  has  evolved  using  somewhat  different  terminology  than  that 
used  up  to  this  point.  The  concepts  in  this  literature  can  also  be  applied  to  other  fields  of 
research.  However,  their  unfamiliarity  to  someone  outside  the  field  can  make  reading  the 
literature  difficult.  The  purpose  of  section  D  is  to  show  how  the  most  important  of  the 
microdosimetric  quantities  are  described  by  the  Two  Component  Model.  At  the  same  time,  the 
rationale  behind  them  will  be  presented  and  they  will  be  related  to  the  more  common  terminology 
in  previous  sections.  D1  discusses  specific  energy  related  quantities  and  D2  discusses  lineal 
energy  related  quantities.  The  specific  energy  and  lineal  energy  are  the  microscopic  analogues 
of  dose  and  LET,  respectively. 


1.  Specific  Energy  Related  Quantities 

a.  Definition  and  Background 

The  term  specific  energy  was  originally  introduced  to  distinguish  between  the  process  of 
dose  deposition  on  a  microscopic  scale  and  the  traditional  concept  of  dose,  which  is  a 
macroscopic,  average  quantity.  The  specific  energy  is  defined  as  the  energy  imparted  by 
ionizing  radiation  to  a  small  volume  divided  by  the  mass  of  the  volume.  This  definition  is 
similar  to  that  of  dose  but  it  is  recognized  that  specific  energy  is  a  quantity  which  refers  to 
localized  dose  deposition  and  has  a  stochastic  nature.  The  word  "localized"  is  emphasized  to 
distinguish  between  the  energy  lost  by  radiation  as  it  interacts  with  a  site  and  the  energy 
deposited  in  the  site.  In  microdosimetry  theory,  specific  energy  is  represented  by  z.  In  previous 
sections  the  symbol  D  has  been  used  to  represent  the  same  quantity  with  an  overlined  D 
indicating  an  average  value.  Specific  energy  distributions  are  often  described  formally  using 
probability  density  functions.  Two  kinds  of  probability  density  functions  have  proven  to  be  of 
interest  in  microdosimetiy.  The  first  weights  a  quantity’s  value  by  the  fraction  of  occurrences 
of  that  value.  This  is  called  a  frequency  weighted  probability  density  function.  The  second 
weights  a  quantity’s  value  by  the  fraction  of  dose  associated  with  that  value.  This  is  called  a 
dose  weighted  probability  density  function.  The  frequency  weighted  probability  density  function 
associated  with  specific  energy  is  denoted  by  f(z)  and  can  be  thought  of  as  an  event  frequency 
spectrum.  The  fraction  of  specific  energy  occurrences  between  z  and  z  +  dz  is  f(z)dz.  The  dose 
weighted  probability  density  function  associated  with  specific  energy  is  denoted  by  d(z)  and  can 
be  thought  of  as  an  energy  deposition  spectrum.  The  fraction  of  dose  deposited  by  specific 
energy  between  z  and  z  +  dz  is  d(z)dz.b  The  reason  for  introducing  the  dose  weighted 
probability  density  function  is  that  the  energy  concentration  in  a  site  appears  to  determine  many 
biological  effects  of  radiation. 


b  In  order  to  avoid  confusion  with  other  notation  involving  the  symbols  f  and  "d", 
probability  density  functions  will  always  be  written  with  their  appropriate  argument. 
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The  more  fundamental  probability  density  functions  are  those  due  to  single  events.  It  is 
generally  recognized  that  because  individual  events  are  statistically  independent,  a  probability 
density  function  for  deposition  of  an  arbitrary  dose  can  be  constructed  from  the  single  event 
probability  density  function.1,2,3,6  The  remainder  of  section  D  deals  with  single  event  related 
quantities  and  their  probability  density  functions. 

Considering  only  single  event  occurrences,  z  represents  the  dose  deposited  in  a  site  by  a 
single  particle.  This  has  been  denoted  in  previous  sections  by  d.  This  is  a  stochastic  quantity 
with  a  frequency  weighted  probability  density  function  fv(z).  The  subscript  "1"  indicates  it  is 
associated  with  single  events.  The  function  fx  (z)  is  relevant  for  analyzing  event  numbers  and 
frequencies  resulting  from  irradiation  of  a  site.  The  fraction  of  single  events  which  occur  in  the 
site  with  a  specific  energy  between  z  and  z  +  dz  is  fx(z)dz.  The  frequency  mean  of  specific 
energy  is  given  by 


ZF 


(119) 


where  the  subscript  F  is  conventionally  used  to  indicate  that  the  average  is  taken  using  the 
frequency  weighted  probability  density.  This  quantity  is  called  the  frequency  mean  specific 
energy  per  event.  It  is  the  average  specific  energy  or  average  dose  deposited  by  a  single  particle 
in  the  site. 

Suppose  a  dose  weighted  distribution  is  of  interest  instead  of  a  frequency  weighted 
distribution.  The  fraction  of  dose  deposited  by  single  events  with  specific  energy  between  z  and 
z  +  dz  is 


dx(z)dz  = 


z-fi(z)dz 


(120) 


The  dose  weighted  probability  density  per  event  is  denoted  by  dx  (z).  It  is  relevant  for  analyzing 
effects  of  local  energy  concentration  in  an  irradiated  medium.  The  dose  mean  of  specific  energy 
is  given  by 
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(121) 


zD  =  fz-dt(z)dz 
o 


The  subscript  D  indicates  that  the  mean  is  taken  using  a  dose  weighted  probability  density.  The 
quantity  is  called  the  dose  mean  specific  energy  per  event.  It  can  also  be  expressed  in  terms 
of fy(z)  by  inserting  (120)  into  (121). 


Jz^'fj(z)dz 

-2 -  (122) 

•a 

fz’f^zydz 


It  can  now  be  seen  that  the  dose  mean  specific  energy  per  event  is  the  ratio  of  the  second 
moment  to  the  first  moment  of  specific  energy  per  event. 


z 


(123) 


Identifying  z  with  d,  this  result  is  equivalent  to  equation  (37). 
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Table  I  summarizes  the  microdosimetry  nomenclature  and  symbols  for  specific  energy 
quantities  per  event,  their  relation  to  notation  used  in  previous  sections  and  their  interpretation. 


TABLE  I.  SPECIFIC  ENERGY  QUANTITIES  PER  EVENT 


Microdosimetry 

Nomenclature: 

Microdosimetry 

Symbol: 

Relation  to  Notation 
in  Prior  Sections: 

Interpretation: 

specific  energy 

z 

d 

dose  deposited  in 
site  by  single 
particle 

frequency  mean 
specific  energy 

ZF 

1 

average  dose 
deposited  in  site  by 
single  particle 

dose  mean  specific 
energy 

Zd 

Ifd 

average  dose 
deposited  in  site  by 
single  particle  using 
dose  weighted 
probability  density 

b.  Specific  Energy  and  the  Two  Component  Model 

In  the  Two  Component  Model,  dose  is  deposited  in  the  site  as  a  result  of  direct  ion  strikes 
and  as  a  result  of  secondary  electrons  created  by  ions  outside  of  the  site.  The  problem  of 
describing  specific  energy  quantities  in  the  context  of  this  model  is  related  to  that  of  a  combined 
or  mixed  probability  distribution.  The  goal  is  to  express  the  overall  frequency  (dose)  mean 
specific  energy  in  terms  of  the  frequency  (dose)  mean  specific  energies  for  ions  and  electrons. 

Recall  that  the  frequency  mean  specific  energy  is  the  average  dose  deposited  in  the  site 
by  a  single  particle.  This  quantity  is  therefore  related  to  the  ion  and  electron  components 
through  equation  (95),  the  expression  for  the  average  dose  of  the  mixed  distribution  of  ion  and 
electron  events.  Rewriting  (95)  in  terms  of  microdosimetry  notation, 


Zf  =  P'ZFJon  +  (1  ~PYzf, 


where  zFion  and  zFe  are  the  frequency  mean  specific  energies  of  ions  and  electrons,  respectively. 
They  are  equivalent  to  the  average  dose  deposited  in  the  site  by  a  single  incident  ion  and  a  single 
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secondary  electron  created  outside  the  site,  respectively.  Substituting  the  expressions  for  P  and 
1  -  P  given  by  (85)  and  (86),  the  following  result  is  obtained. 


N{, 


r=-z. 


Fjton 


N. 


Fjt 


(125) 


The  frequency  mean  specific  energy  can  be  expressed  in  terms  of  an  ion  component  and  an 
electron  component,  which  correspond  to  the  first  and  second  terms  in  (125),  respectively.  The 
ion  component  is  the  average  dose  deposited  by  a  single  incident  ion  weighted  by  the  fraction 
of  ion  events  which  occur  in  the  site.  The  electron  component  is  the  average  dose  deposited  in 
the  site  by  a  single  electron  produced  in  the  surrounding  medium  by  an  ion,  weighted  by  the 
fraction  of  these  electron  events. 

Next,  the  dose  mean  specific  energy  per  event  will  be  written  in  terms  of  the  dose  mean 
specific  energies  for  ions  and  electrons.  Note  that  equations  (37)  and  (123)  both  indicate^  is 
the  ratio  of  the  second  moment  to  the  first  moment  of  the  dose  deposited  by  single  events  in  the 
site.  Equations  (83)  and  (95)  give  the  appropriate  expressions  for  the  second  and  first  moments, 
respectively.  These  were  obtained  from  Deming’s  general  expression  for  the  moment  of  a 
combined  or  mixed  distribution  of  alternative  events  given  by  (70).  Substituting  (83)  and  (95) 
into  (37), 


r  Pdl  +  jl-P)-*!2' 


(126) 


Solving  (96)  for  P  and  (97)  for  1  -  P,  and  substituting  these  expressions  into  the  above  equation 
results  in 


ZD 


(127) 


where  fion  is  the  fraction  of  dose  deposited  in  the  site  by  the  ion  component  and  1  -  fion  is  the 
fraction  of  dose  deposited  by  the  electron  component.  Using  the  expressions  given  by  (100)  and 
(101)  in  the  above,  the  final  result  is  obtained. 
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(128) 


~  f ion  ^DJon 


0  -fuJ  *0, 


This  result  could  also  have  been  obtained  by  comparing  equations  (35)  and  (103).  The  dose 
mean  specific  energy  per  event  can  therefore  be  expressed  in  terms  of  an  ion  component  and  an 
electron  component.  The  ion  component  is  the  dose  mean  specific  energy  for  an  incident  ion 
weighted  by  the  fraction  of  dose  deposited  in  the  site  by  ions.  Similarly,  the  electron  component 
is  the  dose  mean  specific  energy  for  an  electron  weighted  by  the  fraction  of  dose  deposited  in 
the  site  by  electrons.  Expressions  for  the  dose  mean  specific  energy  for  ions  and  electrons  have 
been  discussed  in  detail  and  are  given  by  (104)  and  (105),  respectively. 


2.  Lineal  Energy  Related  Quantities 

a.  Definition  and  Background 

Kellerer  discusses  in  detail  the  limitations  of  the  LET  concept  for  describing  energy 
deposition  by  charged  particles  in  small  sites2,17  One  problem  is  that  transport  of  energy  initially 
in  the  site  to  the  outside,  which  results  from  the  finite  range  of  secondary  electrons,  is 
disregarded.  Another  problem  is  that  LET  is  a  statistical  mean  so  that  energy-loss  straggling  of 
the  incident  ion  is  not  accounted  for.  Other  limitations  are  also  discussed  in  references  2  and 
17.  The  concept  of  lineal  energy  originated  to  distinguish  between  LET,  a  macroscopic,  average 
quantity,  and  the  corresponding  process  which  occurs  in  a  microscopic  site.  The  lineal  energy, 
y,  is  defined6  as  the  energy  imparted  to  a  site  by  one  event  divided  by  the  average  chord  length 
of  the  site,  cl 


(129) 


For  any  convex  volume,  c  is  given  by  4Vq/S0  where  V0  and  S0  are  the  volume  and  surface  area 
of  the  site,  respectively.  The  units  of  y  are  typically  keV/um.  As  was  the  situation  for  specific 
energy,  there  are  two  kinds  of  probability  density  functions  of  interest  in  connection  with  lineal 
energy  -  frequency  weighted  and  dose  weighted. 

The  frequency  weighted  probability  density  function  for  lineal  energy  is  denoted  by  f(y) 
and  it  is  relevant  for  analyzing  lineal  erergy  event  numbers  and  frequencies.  The  subscript  "1" 
is  not  normally  used  here  because  it  is  understood  that  y  is  defined  in  terms  of  a  single  event. 
The  fraction  of  lineal  energy  events  with  magnitude  between  y  and  y  +  dy  is  f(y)dy,  and  the 
corresponding  frequency  mean  of  lineal  energy  is  given  by 
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(130) 


This  quantity  is  called  the  frequency  mean  lineal  energy.  Inserting  (129)  into  (130)  and 
transforming  the  y  variable  to  e  by  noting  that  f(y)dy  =  f(e)de, 


yF 


c 


(131) 


The  quantity  yF  is  similar  to  LET,  assuming  that  the  particle  loses  only  a  small  fraction  of  its 
energy  in  the  site  and  that  escape  of  secondary  electrons  from  the  site  can  be  neglected. 

Suppose  the  effects  of  local  energy  concentration  in  the  medium  are  of  interest.  Then  the 
dose  weighted  distribution  of  lineal  energy,  d(y),  is  used.  Since  y  is  proportional  to  the  dose 
deposited  in  the  site  by  a  single  event,  the  fraction  of  dose  deposited  by  lineal  energy  between 
y  and  y  +  dy  is 


d(y)dy  - 

[y’f(y)dy 


(132) 


The  corresponding  dose  mean  of  lineal  energy  is  given  by 


yD  =  {yd(y)dy 
o 


(133) 


This  quantity  is  called  the  dose  mean  lineal  energy.  Its  relation  to  the  frequency  weighted 
probability  density  function,  f(y),  is  seen  by  substituting  equation  (132)  into  (133). 
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(134) 


/y2  JWy 


fy  'fiy)dy 


Equation  (134)  shows  that  the  dose  mean  lineal  energy  is  the  ratio  of  the  second  moment  to  the 
first  moment  of  the  lineal  energy  distribution  in  the  site. 


yD 


y 

y 


(135) 


It  will  be  recalled  that  equation  (123)  is  a  similar  relation  for  the  dose  mean  specific  energy  per 
event. 


Substituting  (129)  into  (134)  and  transforming  the  integration  variable  to  e,  the  result  can 
also  be  expressed  as 


c  e 


(136) 


48 


Table  II  summarizes  the  microdosimetry  nomenclature  and  symbols  for  lineal  energy 
related  quantities,  their  relation  to  notation  used  in  previous  sections  and  their  interpretation. 


TABLE  H.  LINEAL  ENERGY  QUANTITIES 


Microdosimetry 

Nomenclature: 

Microdosimetry 

Symbol: 

Relation  to  Notation 
in  Prior  Sections: 

Interpretation: 

lineal  energy 

y 

tic 

energy  deposited  in 
site  by  single 
particle  divided  by 
average  chord 
length  of  site 

frequency  mean 
lineal  energy 

57 

tic 

average  energy 
deposited  in  site 
divided  by  average 
chord  length  of  site 

dose  mean  lineal 
energy 

57 

?/(e-f) 

average  lineal 
energy  deposited  in 
site  using  dose 
weighted  probability 
density 

b.  Lineal  Energy  and  the  Two  Component  Model 

First  an  expression  for  the  frequency  mean  lineal  energy  will  be  obtained  in  terms  of  the 
frequency  mean  lineal  energies  for  ions  and  electrons.  Recall  from  equation  (131)  that  yF  is  the 
average  energy  deposited  in  the  site  by  a  single  particle  divided  by  the  average  chord  length  of 
the  site.  In  the  Two  Component  Model,  the  average  energy  deposited  by  an  ion  event  is 


-h 


ion  ^  ion 'S  ion 


(137) 


where  Lion  is  the  ion’s  unrestricted  LET.  Assuming  that  tertiary  electrons  generated  within  the 
site  by  incident  secondary  electrons  are  not  energetic  enough  to  escape  the  site,  the  average 
energy  deposited  by  an  electron  event  is 
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(138) 


where  Le  is  the  unrestricted  LET  of  the  incident  electron. 

From  considerations  similar  to  those  applied  to  equations  (124)  and  (125)  in  the  previous 
section,  equation  (131)  can  be  rewritten  as 


~  .  Eion  + 

F_  * 


C 


(139) 


In  analogy  to  equation  (131),  the  frequency  mean  lineal  energy  for  ions  and  electrons  can  be 
identified  as 


yr. 


JOH 


(140) 


and 


-  — 
c 


(141) 


Then  (139)  can  be  rewritten 


-  N, 


ion 


yF  =  =  _•;» 

*  «  T  ‘  » 


l 


N.  — 


(142) 


Equations  (140),  (141)  and  (142)  are  used  to  evaluate  the  frequency  mean  lineal  energy  in  the 
Two  Component  Model. 

Due  to  the  similarity  of  the  quantities  zp  and  yf,  it  is  worthwhile  at  this  point  to  write 
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down  the  exact  relations  for  these  quantities  and  for  their  components.  Comparing  the 
interpretation  of  TF  given  in  table  I  with  the  expression  for  y}  given  by  equation  (131),  it  is  seen 
that 


1.6x10-“ 

»VPo 


c-yF 


(143) 


In  this  expression,  zF  is  in  rads,  yF  in  keV/pm,  (Tin  pm,  V0  in  cm3  and  p0  in  g/cm3.  Similar 
relations  exist  for  both  ion  events  and  electron  events. 


and 


ZFJon  ~ 


1.6x10  “  - 
- •  c 

K-  Po 


y FJton 


(144) 


1.6xl0'u 

V0  Po 


(145) 


Finally,  an  expression  for  the  dose  mean  lineal  energy  will  be  obtained  in  terms  of  the 
dose  mean  lineal  energies  of  ions  and  electrons.  Comparing  equations  (36)  and  (136),  it  can  be 
seen  that  a  relation  analogous  to  equation  (143)  exists  for  the  dose  weighted  quantities  ~D  and  Yd- 


yD  = 


vo’Po  1  — 
1.6xl0*u  c  Z° 


(146) 


Substituting  (128)  into  the  above  expression  gives  yD  in  terms  of  the  two  components. 


7d--^ 


1 


fion'ZDJan  +  V  ~ fio^ zDa\ 


(147) 


The  dose  mean  lineal  energy  can  therefore  be  written 
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(148) 


‘fion'yDjon  +  V-fiJ-yD* 


where 


and 


VPo 

1.6xlOn 


(149) 


-  V*  Pq  1  - 

^  ■  Tsas^r* 


(150) 


Note  that  these  last  two  relations  can  be  obtained  by  applying  equation  (146)  to  the  appropriate 
type  of  event. 

In  order  to  put  the  last  two  expressions  in  a  form  more  amenable  to  calculations,  they  will 
be  rewritten  in  terms  of  the  appropriate  random  physical  variables.  Inserting  equation  (104)  into 
(149)  gives  the  result 


=  _1  # 

/ 

.4. 

+  F. 

tom 

c 

J  ion 

V 

6^ 

(151) 


Similarly,  substituting  equation  (105)  into  (150)  results  in 


•  =• 
c 


(  I1  «.2  A2 
Le  Se  „  — 

—  •—  +  —  +  F'W. 


(152) 


A  summary  of  the  quantities  on  the  right-hand  side  of  equations  (151)  and  (152)  is  given  in 
section  C2  immediately  after  equation  (105).  Equations  (148),  (151)  and  (152)  are  used  to 
evaluate  the  dose  mean  lineal  energy  in  the  Two  Component  Model. 
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E.  Overview 


This  section  provides  an  overview  of  the  ideas  involved  in  the  Two  Component  Model, 
which  have  been  discussed  in  section  n.  This  will  be  accomplished  with  the  aid  of  flow  charts. 
The  chart  in  figure  4  indicates  the  applied  probability  methods  used  to  combine  the  appropriate 
random  variables  and  variances.  The  sequence  of  operations  begins  with  the  random  variables 
displayed  on  the  left  side  of  the  figure  and  ultimately  leads  to  the  variance  of  the  total  energy 
deposited  in  the  site.  Factors  unimportant  to  the  overview  such  as  conversions  between  energy 
and  dose  and  conversions  between  variances  and  relative  variances  are  not  included  in  the  figure. 

First  of  all,  the  fluctuations  in  energy  deposited  in  the  site  by  a  single  particle  are 
considered.  The  fluctuations  depend  on  the  product  of  the  random  variables  due  to  LET  and  path 
length.  As  shown  on  the  left  hand  side  of  the  chart,  this  applies  separately  to  ion  events  and  to 
electron  events.  Combining  Lion  with  sion  and  Le  with  se  results  in  the  two  indicated  variances. 
(To  simplify  the  notation  in  the  chart,  the  usual  subscripts  "ion"  and  "e”  are  not  included  with 
the  already  subscripted  ”L”  and  ”s”  variables.  Whether  the  variance  refers  to  ions  or  electrons 
is  clear  from  the  argument.  Similar  reasoning  applies  to  the  other  variances  in  the  column  of 
six  variances.)  The  appropriate  expressions  for  these  two  variances  are  obtained  by  applying 
equation  (40)  to  ions  and  electrons,  respectively.  The  key  equation  numbers  are  indicated  at  the 
appropriate  points  in  the  flow  chart.  Next,  energy-loss  straggling  is  considered.  The  variables 
n  and  6  are  combined  making  use  of  Burgess’  variance  theorem  and  that  result  is  modified  using 
the  fact  that  n  is  controlled  by  Poisson  statistics.  Again,  this  is  done  for  both  ions  and  electrons, 
and  the  resultant  variances  for  straggling  are  obtained  by  applying  equation  (44)  to  ion  events 
and  to  electron  events.  If  ionization  fluctuations  rather  than  energy  fluctuations  are  of  interest, 
then  the  random  variables  nIO„  and  r\e  must  be  considered  in  addition  to  the  variables  mentioned 
above.  These  two  random  variables  are  combined  with  their  associated  Fano  factor  to  give  the 
variances  due  to  Fano  fluctuations.  Expressions  for  the  last  two  variances  are  obtained  by 
applying  equation  (50)  to  the  two  situations.  Proceeding  to  the  right  in  the  figure  past  the  first 
six  boxes  which  correspond  to  the  six  operations  discussed  above,  note  there  are  six  resulting 
variances.  Three  contribute  to  the  total  variance  of  a  single  ion  event  while  the  other  three 
contribute  to  the  total  variance  of  a  single  electron  event  in  the  site. 

The  three  variances  associated  with  the  ion  component  are  considered  to  be  the  result  of 
"simultaneous  events"  of  the  type  discussed  in  section  Cl.  They  are  therefore  summed  to  give 
the  total  variance  due  to  single  ion  events,  var(eion).  The  expression  for  this  variance  is  obtained 
from  the  application  of  equation  (51)  to  an  ion  event.  In  a  similar  manner,  the  three  variances 
due  to  electron  events  are  summed  to  give  the  total  variance  due  to  single  electron  events, 
var(te). 

The  next  step  is  to  combine  the  previous  single  event  variances  with  the  corresponding 
random  number  of  events.  This  is  done  with  Burgess’  variance  theorem.  Then,  assuming  the 
event  numbers  are  subject  to  Poisson  statistics,  the  expressions  for  the  variance  of  the  total 
energy  deposited  by  ions,  var(EitJ,  and  the  variance  of  the  total  energy  deposited  by  electrons, 
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var(EJ,  are  obtained.  The  corresponding  expressions  result  from  applying  either  (30)  or  (81) 
to  the  appropriate  event  type.  Up  to  this  point,  the  energy  deposition  by  ions  and  electrons  have 
been  considered  separately. 

The  last  step  is  to  combine  the  variance  of  the  energy  deposition  distribution  due  to  ions 
with  that  due  to  electrons.  The  two  distributions  are  treated  as  "alternative  event"  distributions 
of  the  type  discussed  in  section  Cl .  The  final  result  is  the  variance  of  the  total  energy  deposited 
in  the  site  and  is  given  by  equation  (102). 

A  brief  overview  of  the  construction  of  the  microdosimetric  quantities  zY,  z^,  JY  and  YY 
in  the  context  of  the  Two  Component  Model  will  now  be  given  as  a  conclusion  to  this  section. 
Two  schemes  are  indicated  within  the  flow  chart  displayed  in  figure  5.  The  upper  scheme  shows 
an  overview  of  the  two  frequency  averaged  quantities  while  the  lower  scheme  shows  the  two 
dose  averaged  quantities. 

First,  the  frequency  averaged  quantities  will  be  considered.  To  obtain  the  frequency 
mean  specific  energy  per  event,  zY,  the  first  moments  of  the  quantities  zFion  and  zF  e  are  combined 
using  the  results  of  alternative  event  distributions.  The  result  is  given  by  equation  (125).  This 
is  equivalent  to  calculating  the  average  dose  deposited  by  the  mixed  distribution  of  ion  and 
electron  events. 

Recall  that  specific  energy  and  lineal  energy  are  the  microscopic  analogues  of  dose  and 
LET,  and  are  related  to  each  other  through  equation  (143).  As  indicated  in  figure  5,  the  final 
expression  for  the  frequency  mean  lineal  energy,  %,  can  be  obtained  by  relating  zY  to  yY  and 
incorporating  the  average  chord  length  of  the  target  site.  The  final  result  is  equation  (142). 

Finally,  the  dose  averaged  microdosimetric  quantities  are  considered.  To  obtain  the 
expression  for  zY,  which  contains  the  information  describing  the  single  particle  variability  of 
deposited  dose,  the  combination  of  random  variables  proceeds  as  indicated  in  figure  4  up  to  the 
expressions  for  var(eioJ  and  varfej.  At  this  point  some  algebraic  manipulations  will  convert 
these  two  variances  to  the  expressions  for  the  dose  mean  specific  energy  per  ion  event,  given  by 
(104),  and  the  corresponding  quantity  for  electron  events,  given  by  (105).  These  two  quantities 
are  then  combined  using  the  results  for  alternative  event  distributions  to  obtain  zY,  equation 
(128).  The  expression  for  the  dose  mean  lineal  energy,  ~D,  can  then  be  obtained  by  relating  TD 
to  Yd  which  involves  incorporating  the  average  chord  length  of  the  target  site.  The  final  result 
is  equation  (148). 
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Figure  5.  Overview  of  the  combination  of  the  appropriate  quantities  in  the  Two  Component 
Model  to  obtain  expressions  for  the  indicated  microdosimetric  quantities. 
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m.  SUMMARY 


The  basic  equations  of  the  Two  Component  Model  have  been  derived  in  detail  beginning 
with  well  known  concepts  in  stochastic  theory.  The  model  is  a  first  principles,  analytic 
description  of  energy  deposition  fluctuations  in  ion-irradiated  sites  which  may  be  as  small  as 
nanometers.  The  new  and  key  feature  of  the  model  is  based  on  combining  the  dose  distributions 
which  result  from  two  types  of  single  events  that  occur  in  the  irradiated  site.  The  first  is  a  direct 
ion  interaction  with  the  site.  The  second  is  a  secondary  electron  interaction  with  the  site,  the 
electron  having  been  produced  by  an  ion  in  the  surrounding  medium.  General  expressions  for 
the  frequency  mean  and  dose  mean  specific  energy  per  event  and  for  the  frequency  mean  and 
dose  mean  lineal  energy  have  also  been  obtained  in  terms  of  the  two  components.  All  of  the 
derived  results  are  applicable  to  a  wide  variety  of  conditions,  which  include  different  types  of 
ions,  and  different  target  site  dimensions,  geometries  and  materials.  The  model  can  be  used  to 
calculate  both  the  energy  deposition  and  the  electron-hole  pair  production  fluctuations  in  the  site. 
These  are  all  important  considerations  for  both  experimental  and  theoretical  investigations  of 
small  sites  irradiated  with  ions,  and  has  direct  implications  for  radiation  effects  in  biological 
systems  and  modem  microelectronics. 

The  next  report  will  discuss  how  the  basic  equations  and  quantities  are  evaluated  and  how 
these  results  compare  with  other  published  calculations  based  on  Monte  Carlo  techniques. 
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